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A total of 628 pb~^ of data collected with the ALEPH detector at centre-of-mass energies 
from 189 to 209 GeV is analysed in the search for gauge mediated supersymmetry breaking 
(GMSB) topologies. In this framework, a novel search for six-lepton final states when the stau 
is the next-to-lightest supersymmetric particle (NLSP) and has negligible lifetime is performed. 
Other possible signatures at LEP are studied and the ALEPH final results described for two 
acoplanar photons, non-pointing single photons, acoplanar leptons, large impact parameter 
leptons, detached slepton decay vertices, heavy stable charged sleptons and multi-leptons plus 
missing energy final states. No evidence is found for new phenomena, and lower limits on 
the masses of the relevant supersymmetric particles are derived. A scan of a minimal GMSB 
parameter space is performed and lower hmits are set for the NLSP mass at 54GeV/c^ and 
for the mass scale parameter A at lOTeV/c^, independent of the NLSP hfetime. Including the 
results from the neutral Higgs boson searches, a NLSP mass limit of 77GeV/c^ is obtained 
and values of A up to 16 TeV/c^ are excluded. 
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Outline 



This thesis describes the final searches performed by the ALEPH collaboration at LEP for 
supersymmetry (SUSY) when a light gravitino is the lightest supersymmetric particle. The 
first chapter is devoted to symmetries in particle physics: those of the Standard Model and 
supersymmetry. The limitations of the Standard Model are reviewed and the solutions given 
by supersymmetry are analysed. The fact that SUSY is a broken symmetry, and a very difficult 
one to break, will lead to a variety of possible models. Gravity mediation and gauge mediation, 
the subject of this thesis, will be presented and the production and decay of the predicted new 
particles will be introduced. 

In Chapter ^ the ALEPH detector and its analysis tools are explained. Chapter ^ documents 
the search for selectron or smuon production developed by the author to cover the case of 
heavy neutralinos and light staus with zero lifetime. Chapter ^ summarises all other searches 
carried out by other collaborators within ALEPH in the context of gauge mediated SUSY 
breaking (GMSB) scenarios, including neutral Higgs boson searches. Together with the six- 
lepton topology described in Chapter ^, these results are interpreted in terms of excluded areas 
in the minimal GMSB parameter space. Chapter ^ describes the scan performed by the author 
over these parameters and the combined results on mass limits, cross sections, and the input 
parameters of the theory. 

Chapter ^ will compare ALEPH results with other collaborations and overview the status 
of GMSB models after LEP. A brief look ahead is also given. 



Chapter 1 

The Standard Model and beyond: 
Supersymmetry 

The Standard Model of fundamental interactions is probably one of the most ac- 
curately tested theories in physics. Its foundations will be briefly reviewed here, 
stressing the importance of symmetries as the organising principle behind it. Sym- 
metries provide an important route to advance the understanding and development 
of a theory. But sometimes as well as discovering a fundamental symmetry it is as 
important to know how to break it. The concept of spontaneous symmetry breaking 
is also central to the Standard Model, and it leads to the main untested sector of the 
theory: the Higgs sector. Hereafter, the problems of the Standard Model are outlined 
and a possible solution in terms of a new (ultimate) symmetry will emerge. 
Supersymmetry is regarded as the most likely incarnation of physics beyond the 
Standard Model to be accessible at present experiments. The strong points of the 
theory are studied in this chapter and special attention is given to the consequences 
of the breaking of this symmetry. The different phenomenological consequences for 
current particle colliders will be described paying special attention to the framework 
of this thesis: models where the electroweak and strong interactions (and not gravity) 
mediate the breaking. 
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1.1 Symmetries in physics: the Standard Model 

The Standard Model (SM) of particle physics describes the electromagnetic, weak and strong 
interactions in terms of the internal symmetries they exhibit. The discovery of a symmetry 
in Nature provides the means not only to describe observations, but to construct theories and 
make powerful predictions based on calculations from conserved quantities. 

In the Standard Model [|I], H, |], the forces between the fundamental particles are medi- 
ated by the exchange of spin-1 bosons. Thus the electromagnetic force is 'carried' by the 
photon which couples to charged particles, the weak force couples all fermions by the W and 
Z bosons, and the strong force couples quarks via the gluons and accounts for the stability 
of nucleons. Table lists all four known interactions and the corresponding force carriers. 
The constituents of matter are spin-^ fermions, the leptons and quarks, which come in three 



families (or generations), as listed in Tab. O. 





Carrier particles 


Interaction Strength Range 


Name Mass(GeV/c2) Spin 


Strong as ~ 0.121 lO^^^m 
Electromagnetic ai ~ 1/128 oo 
Weak a2 ~ 10"^ lO^^^m 
Gravitation Gn ~ 10^'^^ oo 


Gluons (g) 1 
Photon (7) 1 
ZO,W± 91.2, 80.4 1 
Graviton (G) 2 



Table 1.1: The four interactions in Nature, with the corresponding boson carriers. The values of 
the coupUngs a are measured at the Z mass; they evolve with energy The graviton has not yet been 
observed. The gravitational constant Gn has units of (GeV/c^)~^. Adapted from Ref. Q. 





Leptons 


Quarks 


Name Q Mass (MeV/c^) 


Name Q Mass (MeV/c^) 


Families 




electron e — 1 0.511 
e-neutrino Uc < 3 x 10^^ 


up u +2/3 1-5 
down d — 1/3 3-9 


2nd 


muon fi — 1 106 
/i-neutrino i^^ < 0.19 


charm c + 2/3 1150-1350 
strange s — 1/3 75-170 


3rd 


tau r - 1 1780 
r-neutrino Ur < 18.2 


top t +2/3 174+5x10^ 
bottom b - 1/3 4000-4400 



Table 1.2: The fermions in the SM are grouped into three famiUes (or generations). Leptons feel the 
electroweak force whilst quarks also feel the strong force. Masses are given in MeV/c^ Only upper 
limits exist on the mass of neutrinos. The quark masses cannot be precisely determined since QCD 
forbids isolated quarks. For each fermion / there is an antifermion / with opposite electric charge 
-Q. 

By imposing certain symmetries, one can write down the Lagrangian and thus derive the 
equations of motion for all particles described above. This is the beauty of the theory: specify 
a global sjTiimetry group and the nature of the interaction is fully determined. So by building 
on experiment, one learns the following symmetries need to be fulfilled by the Standard Model 
Lagrangian: 
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The Standard Model and beyond: Supersymmetry 



• Poincare invariance [^]. The theory must be invariant under continuous spacetime trans- 
formations such as translations, rotations and Lorentz transformations, i.e. the SM 
obeys the laws of Special Relativity. These symmetries imply the conservation of linear 
momentum, angular momentum and 4-momentum, respectively. 

• The CPT symmetry, which swaps the charge, parity and time flow of the process, must be 
respected exactly. This implies the conservation of the corresponding discrete quantum 
numbers. 

• And finally, the gauge (phase) symmetries. The Standard Model is based on the local 
symmetry group SU{3)c ® SU{2)i^ (g) U{1)y- This structure uniquely determines the 
form of the interactions. The SU{3)c symmetry is that of the strong interaction between 
quarks, which can transform one quark to another by means of eight gauge bosons (the 
gluons) conserving the colour charge. Similarly, the gauge invariance under 5'[/(2)l 
requires three gauge fields (eventually the W+,W~ and Z°) and conserves the weak 
isospin II. This symmetry allows an up- type fermion to interact with a down-type 
fermion via the weak force. And the group U{1) is associated with one boson field and 
conserves the hypercharge Y, giving rise to the electromagnetic force. Actually, only 
left-handed fermions which are weak isospin doublets like (e°)L will couple to the SU{2)i^ 
gauge bosons. Right-handed fermions are isospin singlets like ur, and can only interact 
if they carry non-zero hypercharge or colour quantum numbers. 

At this stage the theory is gauge invariant, i.e. transforming the fermion fields in the 
theory by a local phase shift requires the introduction of gauge bosons to leave the Lagrangian 
unchanged. But all the particles in the theory are massless. Indeed, fermion masses should 
appear in terms like mipip = m('?/'L'?/'R -|- ipRiph), which explicitly breaks the SU{2)i^ ® U{1)y 
symmetry in which left-handed fields transform independently of right-handed fields. And 
similarly for vector boson mass terms like ^m^W^W^. Nevertheless, by inspecting Tables |lTl 



and |1.2| , it is obvious that fermions have mass and the weak bosons are (very) massive. The Z 
and W bosons have to be massive to give the weak interactions a very short range. Thus some 
mechanism to break the S'[/(2)l ® f^(l)Y symmetry into t/(l)em is needed if the theory is to 
describe massive W and Z bosons and a massless photon. 

1.1.1 Electroweak symmetry breaking 

An elegant way of introducing massive gauge bosons for the weak interactions without explicitly 
breaking the SU{2)i^ (g) U{1)y gauge invariance of the Lagrangian, was proposed by Higgs and 
others P, 0] using the concept of spontaneous symmetry breaking. Briefly, the introduction of 
a complex spin-0 (scalar) field does the job if it is a doublet in SU{2)i^ and its potential 
respects the most general SU{2)i^ ® U{1)y gauge invariant form: 



where /i is the mixing parameter and A is the quartic scalar self-coupling. For /i^ < and A > 0, 
the above potential has an infinite number of non-zero solutions for which only the field's norm 
is known: |(/)|^ = \{4>i -|- 02 + 03 + ^4) = — /i^/2A = v"^ /2, where v is the energy scale at which 
the electroweak symmetry will be broken. The fundamental state is thus degenerate (the 
circle of vacuum states at V^cp) = in Fig |l.l| b and SU{2)i^ is broken when only one of the 
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hn(^) Re(0) 



Figure 1.1: The scalar potential V{(j)) = ^x^cj)^ cj) + X{(t>^(j))'^ for a single complex field (p = (pi + i(f)2, 
with (a) fj? > and (b) fi'^ < 0. In a the minimum (ground state) is at Kecj) = Imcj) = 0, while for b 
the minimum of the potential is at |(/)| = = 174 GeV, which describes a circle in the (Re(j),lm(j)) 

plane. Any point in the circle is equally likely, thus the spontaneous symmetry breaking occurs when 
only one vacuum state gets chosen, e.g. (pi = v;(f)2 = 0. From Ref. 

infinite points in the circle is chosen by the physical vacuum. Such a choice, for example 
03 = V] 01,2,4 = 0, constitutes a spontaneous breaking of the invariance around the circle, but 
leaves f/(l)em unbroken. This is the essence of spontaneous symmetry breaking in the Higgs 
mechanism: the physical vacuum adopts a specific ground state with non-zero expectation 
value v/\/2 (thus naturally breaking the symmetry) but the Lagrangian remains invariant. 

By spontaneously breaking a continuous symmetry, massless particles called Goldstone 
bosons must be produced 0. In this case the electroweak gauge symmetry had four scalar 
degrees of freedom 0i, 2,3,4, but only SU{2) is broken thus only three massless Goldstone bosons 
are expected. These three degrees of freedom are interpreted as the longitudinal components of 
the three now massive vector bosons W"*", W~ and Z°. The fourth degree of freedom is predicted 
to be the Higgs boson the only particle in the Standard Model yet to be discovered. 

In summary, the Higgs mechanism has allowed the massless weak bosons to become massive 
by eating up the Goldstone bosons arising from the SU{2) breaking, while keeping the photon 
massless. Furthermore, it predicts the existence of a scalar boson with mass mn = -\/— 2/?, 
which is a free parameter and must be experimentally measured. The fermion masses can then 
be generated by Yukawa couplings gf between the fermion fields and the Higgs field of the 
form: g f[ipi^(j)ip^^ + V^r^Vl]- The inner products of Higgs and left-handed fermion doublets 
of SU (2)l yield gauge invariant singlets, so these terms can now be added to the Lagrangian 
without destroying its invariance. 

1.1.2 Defects of the SM 

The Standard Model is not only an elegant description of the symmetries behind matter and 
interactions. It has been tested in detail by experiments and agreement has always been 
found between the measured and predicted values from the theory [||. However, despite its 
experimental success, there remain some theoretical issues which make it impossible for it to 
be a complete description of Nature. They are reviewed here. 
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Gravity 

The Standard Model is an accurate description of the electromagnetic, weak and strong interac- 
tions, but does not include the fourth known force: gravity. Newton's constant has dimensions 
of the inverse of a mass squared (Gn = where Mp ~ lO^^ GeV is the Planck mass where 

quantum effects become important) whereas the other interactions have dimensionless coupling 
constants. In General Relativity, gravity, as a force, manifests itself through deformations in 
spacetime. In the Standard Model the spacetime is assumed to be flat and forces to arise from 
the exchange of quanta. This is probably one of the remaining great issues to be resolved in 
physics: how to make the Standard Model and General Relativity compatible. 



Number of free parameters 

Ideally, one would expect that the three coupling constants (ai, 0^2, as for each symmetry group 
in Tab. |1.1| ) would be the only free parameters of the theory. Assigning the matter particles 
to a specific representation of the gauge group and then measuring the couplings should be 
enough to have a complete theory of Nature. Conversely, the SM has at least 28 arbitrary 
parameters to account for: three gauge couplings; two parameters in the Higgs potential (/^ 
and A); 12 fermion masses (including the neutrinos); three mixing angles and one phase in the 
quark sector and another three mixing angles and three phases in the lepton sector for massive 
neutrinos; and finally one CP violation phase in the strong interaction. 



Fermion masses 



Looking back at Tab. |1.2| the diversity of observed masses for the fermions is striking. As was 
mentioned in the previous section, the Higgs mechanism incorporates massive fermions into the 
theory while keeping the Lagrangian gauge invariant. The couplings between the Higgs and the 
fermion fields are proportional to the fermion masses: gf = V2mf/v. If v/V2 = 174 GeV, why 
is the top-quark Yukawa coupling so close to one and the electron's five orders of magnitude 
weaker? The SM does not explain this mass hierarchy. 



The recent experimental evidence for neutrino oscillations |]T0|, |TT|, only possible if they 
are massive, strongly suggests physics beyond the SM. The extreme smallness of their masses 
{rriiy < 0.2eV/c^ ~ 10~^me) indicates that their origin cannot be the Higgs mechanism. Such 
small masses may only arise if the existence of right-handed neutrinos is postulated (Dirac 
masses, e.g. VktyidVi^) or if a Majorana nature of the neutrinos is invoked (if the neutrino is 
its own antiparticle: vj^mvi,) [|l^. In both cases, the origin of neutrino masses must be linked 
with an energy high enough to violate lepton number L. 



Dark matter and the baryon asymmetry 



Conventional matter, made of baryons like the proton and the neutron 
of the total energy content of the Universe according to Big-Bang nucleosynthesis 



represents < 5% 
r3[ and 



cosmic microwave background (CMB) radiation |14| analyses. Studies of large-scale structure, 
rotation curves in spiral galaxies, luminous matter and CMB, suggest that the total amount 
of dark matter in the Universe corresponds to about 30% of the total energy density pc ~ 
5 X 10-6 GeV/ cm3 



15 



If the SM seems to describe only 5% of pc, what kind of matter is this 
remaining 30%? 

Furthermore, why is the Universe made of matter and not antimatter? The SM equations 
are the same for matter and antimatter, but then there is no evidence for antimatter galaxies 
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or stars in the Universe, only matter. Three conditions are required to generate a baryon 



asymmetry |T6|: the baryonic quantum number B cannot be conserved; the C and CP sym- 
metries are violated; and a period of thermal inequilibrium is necessary for the asymmetry 
not to be washed away. Although the SM could accommodate in non-perturbative processes a 
violation of the baryon number, and CP violation arises naturally in a three generation model, 
it seems the amount of CP violation predicted in the SM is not sufficient [|I7|. Today there 
are ~10^ photons for each baryon in the Universe, so at some point there had to be roughly 
one extra quark per billion qq pairs. The origin of this asymmetry requires the existence of 
physics beyond the SM with a stronger CP violation. 



Hierarchy or naturalness problem 

Paradoxically, it is a particle whose existence is not yet proven that poses the most serious 
threat to the Standard Model. The Higgs boson mass is a free parameter in the Standard 
Model, but for the theory to be valid it must be bound from below and, specially, from above. 
However, the Higgs mass in the Standard Model gets contributions from loops with the gauge 
bosons, fermions and with itself, as seen in Fig. |1.2|. The problem arises when considering 



radiative corrections from scalar particles, as in Fig. |1.2| c. Gauge boson and fermion masses 
are protected against divergences, so their loops do not pose a problem. 





(a) W and Z boson contribution 



(b) Fermion contribution 



(c) Scalar contribution 



Figure 1.2: One-loop quantum corrections to the Higgs mass. 



Scalar particles tend to get masses of the order of the largest mass scale in the theory. 
And if the SM is to hold unscathed up to the energy where gravity becomes strong, Mp, then 
the quantum corrections to the Higgs mass diverge. Normally one introduces a cutoff energy 
scale A between the W mass and Mp, above which the SM is invalid and some new physics 
must appear. Thus the Higgs mass will be given by some 'bare' mass term mg (where no 
assumption on new physics is made) and the corrections that involve the new energy scale: 
~ itIq + drrijj- -|- . . ., with: (5m|^ = (9 A^. Then, if on the one hand the Higgs mass 
involves radiative terms of order A^ ~ 10^^ GeV^, but on the other hand it is known to be of the 
order of the electroweak symmetry breaking scale (m^ ~ —v'^), how is it possible that the bare 
mass is so finely tuned as to exactly absorb the enormous radiative contributions and be left 
with the correct physical value? This unnatural fine tuning (of up to 34 orders of magnitude) 
is referred to as the hierarchy or naturalness problem. It can also be expressed as: why is the 
W mass much smaller than the Planck mass? Or equivalently, why is the Coulomb potential 
for a typical particle mass m in the atom so much greater than the Newtonian potential: 
e^ > Gn^^ = mVMp2? 

All these reasons lead to the idea that the SM is only an effective theory at the present 
energy scale ~100 GeV and that there must be a new more fundamental theory just around the 
corner, at A < 0{1) TeV. Supersymmetry is widely accepted as the most promising extension 
of the SM to explain some of the theoretical problems outlined here. 
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1.2 The last symmetry: Supersymmetry 

The symmetries of the SM involve scalar charges, such as colour or hypercharge, that link 
particles with the same spin: Q\J) = \J). The basic idea in supersymmetry (SUSY) is to have 
fermionic charges that relate fermions and bosons, so that: 

Qlfermion) = |boson) and Q|boson) = |fermion). (1-2) 

Supersymmetry is therefore a transformation that relates states of different spin. Particles like 
bosons and fermions that in principle have a completely different nature may be linked and 
form part of the same representation. 

Perhaps the most compelling argument for supersymmetry is that it is the last undiscov- 
ered symmetry allowed in scattering processes. Considering possible extensions of the Poincare 
group, Coleman and Mandula |TB[ proved that the addition of any new operator which trans- 
forms as a boson leads to a trivial S-matrix, i.e. in particle scattering experiments the only 
allowed outcomes would be completely forward or backward scattering, which of course com- 
pletely disagrees with observation. But this very strict 'no-go' theorem can be avoided if the 
operator is fermionic and as was demonstrated later |jl9|, supersymmetry is actually the only 
extension of the Poincare group which does not lead to a trivial S-matrix. So this argument 
should be persuasive enough: 'in Nature something not illegal is compulsory'. 

Supersymmetry is not only attractive aesthetically as an underlying theory for the SM, 
it is also a step towards unification with gravityQ Furthermore, it cures several of the SM 
shortcomings. 

The supersymmetric algebra derived from the above conditions requires the existence of 
new particles with the same mass, gauge quantum numbers and couplings as the SM ones. 
Thus for each fermion in the SM there is a bosonic partner and vice versa, structured in 
supermuUiplets like: 

/ fermion( J = |) ^ f &■ boson( J = 1) \ f graviton( J = 2) \ q o\ 

\^ sfermion( J = 0) J ' \^ gaugino( J = |) J ' Ygravitino(J = |) J 

which describe matter and Higgses, gauge fields and gravity, respectively. The convention in 
naming the superpartners (or sparticles) is to add a prefix 's' (for scalar) to each fermion and 
a suffix 'ino' for each boson. Could any of the known particles be one of these superpartners? 
No, it is impossible to pair together any of the known fermions with the known bosons, for the 
simple reason that their internal quantum numbers do not match. For example, leptons have 
non-zero lepton number L = 1, but bosons have L = 0, thus they cannot form part of the same 
supermultiplet. 

By predicting this plethora of new particles supersymmetry solves neatly the hierarchy 
problem. The radiative corrections from fermions and scalars (diagrams b and c in Fig. |1.2| ) 
now cancel each other naturally at all levels of loops . For each positive scalar loop there is 



now a negative fermion loop with the same couplings. The radiative contribution to the Higgs 



mass becomes 20 : 



,2 



^Making supersymmetry a local (spacetime dependent) symmetry necessarily involves gravity. Local SUSY 
transformations are equivalent to a local coordinate transformation; therefore supersymmetry is entangled with 
a new approach to spacetime, which must be present in any viable 'Theory of Everything' pO|. 
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~ O {ml - m)) (1.4) 

where the quadratic divergences have disappeared and the residual term is much smaller than 
w?jj if the masses of the SUSY bosons and fermions are similar: 



-m5|<lTeV2 (1.5) 



This means that a light Higgs boson, such as electroweak data constraints I^Sj and recent direct 
searches seem to suggest ||2^, might be only possible if SUSY exists at a relatively low energy. 

Another appealing feature of SUSY is that it permits the unification of the three forces at a 
very high energy scale, which cannot be achieved within the SM. The three coupling constants 
(«!, q;2, cts) are not constants at all, but depend on the energy. Every particle with mass nix 
below a scale Q will contribute to the running of the couplings via loop corrections to the force 



carriers propagators and vertices p4| : 



^ ^ +87r6iln(— ) i = l,2,3 (1.6) 



ai(Q) ai{mx) \mx 

where the coefficients 6j depend on the number of colours and the number of active flavours 
(particles whose mass threshold is below Q), including Higgs bosons. These coefficients are 
different in the SM and SUSY, as can be seen in Fig. |1.3[ In the SM the three lines do not 
intersect at a same energy, even allowing for experimental error bands on the measurements at 
present energies. But SUSY models do in fact permit a complete unification of the coupUngs 
at around Qu = 10^^ GeV. Although this is a nice feature of supersymmetry, it does not take 
into account anything else happening during some 14 orders of magnitude from the scale where 
SUSY is turned on up to the GUT (Grand Unified Theory) scale. 

There are two more characteristics of SUSY models worth discussing. The first is that two 
Higgs doublets are needed to give mass to up-type fermions and down-type fermions separately. 
In the SM only one doublet and its complex conjugate (p^ was sufficient, but in SUSY models 
terms are not allowed since they do not transform appropriately under supersymmetric 
transformations. Then, two Higgs doublets 0i = {4>i,4>i) 02 = {4>t:'p2) with opposite 
hypercharge (=Fl) are needed in order to give masses to the SM fermions. Both will acquire a 
non-zero vacuum expectation value, fixed like before: f ^ + f | = (246 GeV)^, but their ratio is 
a free parameter, called tan/?: 

, o _ (02) V2 , 

tan^ = — — = — 1.7 

(0l) 

This makes phenomenology somewhat more interesting since there are now eight degrees of 
freedom and only three of them get 'eaten' to give mass to the weak gauge bosons, leaving 
five physical massive particles: the lighter CP-even neutral h, the heavier CP-even neutral H, 
the CP-odd neutral A and two charged H^. Once the two Higgs supermultiplets have been 
introduced, they can couple to each other in terms like /i0i02; where /i is the mass mixing 
parameter, analogous to the Higgs mass ttih in the SM. 

The other very important phenomenological consequence of SUSY models is that new 
terms which violate lepton- and baryon-number are now allowed in the Lagrangian and still 
are consistent with gauge symmetries. Couplings which violate both B and L are strongly 
constrained from nuclear physics and specially from the proton lifetime [^. In the SM, B— L 



conservation was never imposed by hand but rather was an 'accidental' outcome of the theory. 
Since the strength of these couplings is severely limited, one can impose the conservation of a 
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Figure 1.3: Evolution of the gauge couplings Oi as a function of the energy Q, with the particle 
content of the Standard Model (above) and with the minimal supersymmetric extension of the SM 
(below). The energy scale where SUSY appears has been assumed to be 10^ GeV. 



new multiplicative quantum number: R-parity, defined as: 

Rp = (_1)3B-L+2S ^^_g) 

where S is the spin of the particle. Standard Model particles will have Rp = +1 and SUSY 
particles Rp = —1. To assume R-parity conservation, which is to some extent to impose B— L 
conservation by hand, has three fundamental consequences: 

• sparticles are always produced in pairs, 

• heavier sparticles decay to lighter ones, and 

• the lightest sparticle is stable because it has no allowed decay mode. 

The last point is perhaps the most attractive, since it implies the existence of heavy stable 
particles. If these are neutral they could very well constitute that 30% of the energy density of 



the Universe that must be in the form of cold dark matter ||T^ . However, R-parity conservation 
is not obligatory and there may exist terms, with maybe significant strength, that only violate B 
or L. Nonetheless, the work throughout this thesis assumes that R-parity is strictly conserved. 



1.2 The last symmetry: Supersymmetry 
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g 


gluino 
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1 
1 




7 
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zino 

neutral Higgsinos 


\ 

2 

X 

2 
1 
2 




Al,2,3,4 

neutralinos 


W± 


1 




wino 


1 
2 


Xt2 









charged Higgsinos 


1 
2 


charginos 


G 


2 


G 


gravitino 


3 
2 


G 



Table 1.3: The MSSM particle zoo. 



1.2.1 Supersymmetry breaking 

The major problem with SUSY is that it must be a broken symmetry. There is not a single piece 
of experimental evidence for SUSY in particle physics so far. If a supersymmetric companion 
of the electron existed with its same mass and spin-0 it would have been discovered long ago. 
The problem then is to write a SUSY Lagrangian with terms that make sparticles heavier 
than their SM counterparts, but do not introduce quadratic divergences to the Higgs mass 
(Eq. p..5[). These terms are thus said to induce a 'soft SUSY breaking'. The origin of such 



terms is not known and several possibilities exist So in general one writes all possible 

terms and does not make any assumption about how exactly SUSY is broken. But to 
parametrise our ignorance on the origin of these terms, at least 105 new free parameters 
must be introduced. They account for the sparticle masses, their mixing angles, 40 new CP 
violating phases, etc. . . Thus a theory that was simple and extraordinarily predictive becomes 
almost unmanageable, plagued with many more free parameters than before. The hope is 
that these parameters will be eventually explained by GUT-scale physics in terms of a few 
fundamental ones. Making some assumptions on how SUSY breaking is transmitted and 
imposing experimental constraints will lead to more predictive models as will be seen shortly. 

The minimal supersymmetric extension of the Standard Model, or MSSM for short, can 
then be constructed in the most general way, allowing for R-parity violating terms and soft- 
SUSY breaking terms. The properties of the mass spectrum are now reviewed. 

Table |1.3| lists the complete particle content of the MSSM. For each left- or right- handed 
fermion /l,r there exist two different sfermions /l,r with spin-0 and weak isospin ^ and 0, 
respectively. Of course the sfermions being scalars do not have a defined helicity, but they do 
have different weak isospin depending on the helicity of their SM partners. Once electroweak 
symmetry is broken these different weak isospin states are allowed to mix. Thus for each 
flavour a 2 X 2 mixing matrix exists which will give the physical mass states from the weak 



eigenstates /l,r. It takes the following general form pO| : 

i2 I ^2 I ^2 



/l /e) ' t ,„2 , L , ,„2 \ (1-9) 



'r 
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with 

ml) = m| cos 2(3 (Js + Q sin^ 6w) 

tan/3\ r f e,iM,T,d,s,b (1.10) 



= (^/ - /^cot ^) for / 



u,c,t 



where I3 is the third component of weak isospin, sin^ 9w = 0.23 and Af are the trihnear Higgs- 
sfermion-sfermion couplings arising from soft SUSY breaking. The off-diagonal terms depend 
on ruf, so it is clear that there will be more mixing for ^lr, and if tan/? is large also for 6l,r 
and tl r. The mass eigenstates can be obtained by diagonalisation. 



fi \ ^ f cos 6 J sin 6 J \ f fh 
f2 J \ sin 6 J cos Oj^ J \fR 



[1.11) 



are 



where Oj is the mixing angle. For example, in the case of stau mixing, Eqs. |L9| and |1.11 
then related by: 

rrir tan 29 f 

In general, the lighter mass eigenstate fi will correspond to the right handed tr state if there 
is no mixing {Of = and = fi tan (3). As the mixing increases (i.e. for increasing values 
of the mixing angle) the lighter fi becomes an admixture of right- and left-handed staus, 
until 6f = 45° which corresponds to maximal mixing. In this case the right- and left-handed 
components of the lightest stau are equal. Mixing is therefore more relevant in the third family: 
the lightest stau can be lighter than selectrons and smuons, and the lightest stop will generally 
be the lightest squark. Mixing will also occur between the massive fermionic companions of 
the electroweak gauge bosons. The photino 7, zino Z and neutral higgsinos Hi 2 will mix to 
form the 'neutralinos': Xi 2 3 4; where the lower the index the lighter the particle. Similarly, the 
charged wino and higgsino give four charged mass eigenstates: the 'charginos' Xi2- 
The corresponding mixing matrices for neutralinos and charginos can be found in Ref. p 



1.2.2 Super-Higgs mechanism and the gravitino 

Returning to the problem of SUSY breaking, it turns out to be impossible to break SUSY 
in a phenomenological acceptable way if the only particles and interactions are those of the 
MSSM. In the SM, the Higgs vacuum expectation value determines the scale of electroweak 
breaking. But the specific masses of the bosons and fermions are dictated by the coupling 
of the forces that communicate the information of electroweak breaking: gauge and Yukawa 
couplings, respectively. In the MSSM, to give masses to the gauginos, for example, soft terms 
in the Lagrangian should contain tree level interactions of the type scalar-gaugino-gaugino, but 



these are not allowed in supersymmetry |2^. Similarly, squarks would be unacceptably light 
if their mass was generated at tree level as is the case in the SM . 

Thus the soft mass terms must arise indirectly or radiatively: from the coupling of new 
particles to the particles in the MSSM. Supersymmetry breaking must therefore occur in a 
'hidden' sector at a large energy scale. There the spontaneous breaking of local supersymmetry 
occurs and as a result some field F condensates, acquiring a non-zero expectation value Fq. 
Similarly to the Higgs mechanism, a Goldstone degree of freedom must now appear, but in this 
case it is a fermion since SUSY is a fermionic symmetry. This massless particle is the so-called 
'goldstino', and its two polarisation states are 'eaten' by the massless gravitino, giving it a 
total of four polarisation states to become a massive spin-| particle. The end result is then 
a massive gravitino and a massless graviton: local supersymmetry is thus successfully broken. 



1.3 Gravity mediated SUSY breaking 
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This super-Higgs mechanism is actually the only consistent way of breaking SUSY, as the 



Higgs mechanism was the only consistent way of breaking gauge symmetry . The gravitino 
mass is then given by: 

rrir = — (1-13) 
^ V3Mp 

where Fq is the energy scale at which SUSY breaking occurs in the hidden sector (like f ^ in the 
Higgs mechanism) and Mp is the reduced Planck mass: Mp = (SttGn)"^^^ = 2.4 x 10^^ GeV/c^. 
By dimensional analysis, it is intuitively clear that the gravitino should be massless in unbroken 
SUSY (Fo 0) or if gravity becomes irrelevant (Mp — oo or Gn — ^ 0). 

Now MSSM particles still have to be informed of the SUSY breaking. This can be achieved 
in different ways. Since all massive particles feel gravity, a natural candidate to communicate 
the hidden sector with the 'visible' sector where all MSSM particles sit, is the gravitational 
force. Thus gravity mediated SUSY breaking or supergravity (SUGRA) models, postulate 
the coupling between Planck-scale physics and the MSSM particles via gravity (exchanging 
gravitons). These weak couplings (of the order 1/Mp) generate acceptable soft mass terms and 
a definite MSSM spectrum. Another possibility is that the SUSY breaking is transmitted to 
the MSSM particles by the usual gauge SU{3)c ® SU{2)i^ ® U{1)y interactions. This case is 
referred to as gauge mediated SUSY breaking (GMSB) models and they imply the existence 
of an intermediate energy sector of particles, called 'messengers' $ which will directly feel the 
SUSY breaking and then transmit it to the MSSM particles by radiative corrections, giving 
rise to a given MSSM mass spectrum. This last type of mediation is the framework of this 
thesis. A schematic drawing of the structure of SUSY breaking and mediation for these two 
scenarios is presented in Fig. WJ. Both alternatives are discussed in the next sections. 



A very nice feature of the super-Higgs mechanism is that it naturally induces electroweak 
symmetry breaking. If one assumes, as the unification of couplings suggests, that all fermion 
and boson masses unify at the GUT-scale with values mo and mi/2 respectively, one can 
calculate all MSSM masses in terms of these GUT-scale input values. It turns out that mg 
evolves with energy in such a way that drives the squared mass of the up-type Higgs negative 
due to the large (negative) radiative corrections from the top-quark Yukawa coupling and the 
stop mass. Indeed, electroweak symmetry breaking is generated radiatively and a link between 



the SUSY breaking Higgs masses and the electroweak scale mz can be obtained ||29 



^ ~ 2 + tan2/3-l ^'''^^ 

This again expresses the hierarchy problem in a different format: the n parameter on the left 
can be as big as the highest energy in the theory (it does not break any SM symmetry), but 
both mz and the Higgs masses on the right have to be of the order the electroweak scale. What 
makes n then lie in a physically acceptable region? 



1.3 Gravity mediated SUSY breaking 

If gravitational interactions are responsible for the soft mass terms, Planck-scale physics is 
the messenger of SUSY breaking . A simplifying assumption, although not inevitable, is to 
consider the universality of masses and couplings. First, since the gauge couplings seem to unify 
at Qu = 10^^ GeV (see Fig. pT^ ), the gaugino masses Mi, M2 and M3 for the photino, wino and 
gluinos respectively are also assumed to unify with a common value of: mi/2 = Mj ~ tuq [pO|] . 
Secondly, the scalar masses-squared and the trilinear couplings are taken to be flavour- diagonal 
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Figure 1.4: Supersymmetry breaking and its mediation. SUSY is broken (spontaneously) in a high- 
energy 'hidden' sector and its effects are transmitted radiatively to the 'visible' sector where all MSSM 
particles reside. The mediation can occur via gravity or gauge interactions. 



and universal at ~ and Af = A^Yf ~ ttLq, where Yy are the Higgs-fermion Yukawa 
coupling matrices. All the above relations are evaluated at the unification scale Qu or Mp, 
and all of them are of the order Fg/Mp set by gravitational interactions. To calculate the 
mass of each particle in the visible sector, one runs the energy down from the known energy 
dependence of these parameters p8[ . This set of assumptions lead to a very predictive model, 
called minimal supergravity (or mSUGRA), where only five parameters are needed to fully 
calculate the mass spectrum in the MSSM: 

mi/2, mo, Ao, tan/3, sign(/i) (1.15) 



where the absolute value of // can be calculated with Eq. |1.14| and only its sign remains a free 
parameter. 

Since the soft masses have to be of the order of lOOGeV/c^, and these are determined by 
parameters proportional to Fg/Mp, the scale of SUSY breaking a/Fq in these type of models is 
predicted to be of the order of 10^^ GeV. In this type of models the gravitino is then expected 
to be heavy, comparable to the MSSM particles with ~ lOOGeV/c^ or so. 



It is then found that the gaugino masses at any scale satisfy the relation |28 



Q 

M3 = — sin^ ewM2 = cos^ OwM^ (1.16) 
a oa 

which at the electroweak scale (0:3 = 0.118 and a = 1/128) becomes: M3 : M2 : Mi ^7:2:1. 



1.4 Gauge mediated SUSY breaking 
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It can also be derived that in mSUGRA the neutralino Xi is usually the lightest supersymmetric 
particle (LSP) [^. If R-parity is conserved, signatures for SUSY in colliders will consist of 
missing energy carried away by the heavy, stable and very weakly interacting neutralino. 



1.4 Gauge mediated SUSY breaking 

In GMSB models the splitting of masses in the MSSM is generated not at the Planck scale, 
but at some lower energy scale where pairs of very heavy messenger quarks and leptons exist. 
The gravitational interaction is still responsible for the communication between the hidden 
sector and this messenger sector, but its effects are now negligible in the MSSM. Gravity will 
couple to a singlet S in the messenger sector which will acquire a vacuum expectation value 
\/Fin. This singlet breaks supersymmetry in the messenger sector and gives mass to all other 
messenger superfields $ there. In turn, these messengers are able to communicate with the 
MSSM via gauge interactions, but now from a much lower scale than Mp. 

The key element in these models is then the structure of the messenger sector. The heavy 
messengers can be generally described as supermultiplets of SU{5) D SU{3) <S) SU{2) ® U{1), 
for which the parameter N5 is introduced and represents the number of families of To a 
good approximation, once SUSY is broken in the messenger sector, the fields $ will share a 
common mass Mm. 

The MSSM masses are then generated at one loop in the case of gauginos Aj and two loops 
in the case of scalars /, taking into account that couplings are in both cases given by gauge Oj 
strengths. Examples of the Feynman diagrams contributing to the MSSM masses are shown 
in Fig. 




(a) Gauginos (b) Scalars 

Figure 1.5: In GMSB, the (fermionic and bosonic) messengers ^ give supersymmetry breaking 
masses to (a) gauginos A at one loop and to (b) sfermions / at two loops. The coupling is via gauge 
interactions. 



The masses can then be calculated at the messenger mass scale M^, as p9| : 

(1.17) 

A^Ns (1.18) 

where the parameter A = Fm/Mm represents the universal mass scale of SUSY particles (like 
mo or mi 12 in SUGRA). The most noticeable consequence is the degeneracy between squark 
and slepton masses. In SUGRA models the different flavour masses arise directly from a 
high energy Mp where flavour symmetry is respected. There is no control then on the mass 
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difference between squarks and sleptons, which could be in principle very large. These possible 
flavour-breaking soft mass terms are very dangerous, because they lead to flavour changing 
neutral currents (FCNC) and large CP violation. Both of which are strongly constrained by 
observations [^]. Measurements of the K^-K^ mass difference and — > e7 decays set stringent 
bounds on the possible mass difference between squarks and sleptons In GMSB models 



instead, sfermion masses arise at a lower energy A, much below the scale at which flavour 
symmetry is expected to be valid. Thus the SM Yukawa couplings are already present when 
the mass generation takes place, and FCNC are then automatically suppressed. Also, the 
trilinear couplings Aj are negligible at the scale, since they are generated at two loops and 
are further reduced by a factor aijAm. They will be in any case very small at the electroweak 
scale and are not considered as a relevant parameter. 

Another important characteristic of these models is the gravitino mass. Since soft masses 
are proportional to the coupling constants, the LSP will be the gravitino, the least interacting 
sparticle. By reducing the energy scale at which soft mass terms are generated from Mp to 
Mm, the 'effective' SUSY breaking scale (in the messenger sector) \fW^ is now much lower. 
The messenger mass could be ~ lOOTeV/c^, thus to produce MSSM particles with the right 
mass (m ~ Fm/Mm ~ lOOGeV/c^), v^F^ can be as low as lO^GeV/c^. From this 'effective' 
SUSY breaking scale the gravitino mass is derived now as p9| : 

2,4 ' 



^ wis; ^ T (Mifv j ""1^ (^■^^) 

where Fm = fcFo and is a model dependent parameter describing how SUSY breaking is 
communicated from the hidden sector to the messenger sector. 

If the gravitino is the LSP all MSSM particles will eventually decay into it, and since it 
couples gravitationally one would naively expect those decays to be extremely slow. But the 
gravitino contains two longitudinal components from the goldstino it absorbed by means of 
the super-Higgs mechanism. And these components with helicity ±^ and gauge couplings may 
lead to decay rates high enough to be of experimental importance. 

Finally, to completely specify the MSSM spectrum and its phenomenology, it is only nec- 
essary to use six parameters: Mm, m^. A, N5, tan/5 and sign(/i). Hence, GMSB models are 
highly predictive as will be discussed in Chapter |[ 



1.5 Collider signatures 

The best probe for supersymmetry is naturally direct production of superpartners at high 



energy colliders. Figure |1.6| shows the Feynman diagrams for the production of all relevant 
SUSY particles at LEP, where electrons and positrons collide head-on. 

In Gauge Mediated SUSY Breaking models, experimental signatures differ signiflcantly 
from models with a neutralino LSP. If R-parity is conserved, as is assumed throughout this 
work, all SUSY decay reactions will terminate in the next-to-lightest supersymmetric particle 
(or NLSP) decaying to its SM partner and the gravitino LSP. The identity of the NLSP, 
therefore, is central to the phenomenology. Furthermore, the NLSP may be long lived and 
thus very striking signatures are possible. The NLSP decay length is controlled by the SUSY 
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(b) Neutralino and chargino production 




(d) Higgsstrahlung and associated production 



Figure 1.6: Production of SUSY particles and MSSM Higgs neutral bosons at LEP when R-parity 
is conserved. 
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(a) Contours of the NSLP decay length (Eq. |l.2C 
with = 1) as a function of the NLSP mass and 
the gravitino mass, for pair-produced particles at 
E = 104 GeV. 
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NLSP as a function of the number of messenger 
families and the messenger mass. The solid lines 
represent different values of tan /3 = 3, 10, 30. From 
Ref. (31 



Figure 1.7: The NLSP nature as a function of parameters. 



breaking scale or equivalently (Eq. |1.13|) by the gravitino mass: 

0.01 /100GeV/c2^ 




Anlsp = CT^LSPlP = -^y ^j^^gp J ( o .„w/„2 ) \l—2 ^ (1.20) 

where k,^ is the photino component of the neutralino NLSP, or one in all other cases (see 
Fig. |1.7] a). If rriQ is large enough and the NLSP decay length is measured as a vertex displace- 
ment, it provides a unique method to probe the value of v^Fq, the scale of SUSY breaking, as 
opposed to other measurements which in general are only sensitive to A. Searches for NLSP's 
decaying in the middle of the detector also benefit from lower background rates, hence the 
interest in studying these signatures. 

As regards the NLSP type, depending on the GMSB parameters it will generally be either 
the neutralino or the lightest stau (see Fig. p77[ b. This can be deduced from Eqs. |1.17| and |1.18| 



where sparticles with only U{1)y interactions will be lighter than the rest since ai < ^2 ^ c^a- 
The tiny selectron and smuon mixing factors (Eq. |1.10D allows one to generally treat them as 



degenerate, unmixed mass eigenstates, thus mg = will be assumed for the rest of this work. 
In some cases, specially for low tan (3, the three right handed sleptons are degenerate in mass 
and act as co-NLSP. Also possible in a much more restricted corner of parameter space is the 
neutralino-stau co-NLSP where the mass difference between the neutralino and the sleptons 
is less than the tau mass. Table |1.4| lists the conditions in the mass hierarchy for each type 
of NLSP scenario and its decay |^2|]. A neutralino NLSP will predominantly decay to 7G and 



not ZG for < lOOGeV/c P3[. From now on, references to the slepton NLSP case or 
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Table 1.4: The four possible NLSP scenarios, with the corresponding NLSP decay, neglecting the 
electron and muon masses. 



scenario, should be understood as corresponding to the stau NLSP region, including also the 
more restricted case of sleptons as co-NLSP. 

The superpartner mass spectrum in GMSB models will ultimately determine which SUSY 
particles can be produced at LEP and their subsequent decay chains. The strongly interacting 
states, including the generally hghter stops, are usually too heavy to be produced at LEP. 
Only the lightest neutralino x?? slepton £^ or Higgs will in general be accessible at LEP2. 
The charginos, heavier neutralinos and left-handed sleptons are too heavy to be relevant to 
discovery at these energie^. A typical mass spectrum of GMSB models is shown in Figure |0 



In the following, the lightest neutralino x5 will be abbreviated to x ci^nd referred to as 'the' 
neutralino, when no other heavier neutralino or chargino is present. The symbol i will stand 
for the three charged leptons, e, /i or r; and / will only refer to the electron or muon. 




Messenger scale (GeV) 



Figure 1.8: Sparticle masses dependence with Mm in a GMSB model with N5 = 1, tan/? = 3 and a 
messenger scale bino mass of 115GeV/c^. From Ref. |]3^ . 
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Even if no other SUSY particle is accessible at LEP, direct production of gravitinos is still possible. The 
cross section for e+e" GG7 scales with 1/m^, so only ultra light gravitinos could be observed through 
direct production. ALEPH has searched for this process and derives an upper limit on the gravitino mass of 
1.3 X lO-^eV/c^ @. 



Chapter 2 

Experimental Overview 



The Large Electron Positron collider (LEP) started operation in 1989 and ended in 
the fall of 2000. The largest accelerator of its kind in the world, it was built to 
study in detail the Z and W massive vector bosons. Prom 1989 to 1995 it ran with 
a centre-of-mass energy (s/s) close to the Z mass, scanning the electroweak sector 
and the properties of the Z particle. This phase is known as LEPl. Prom 1996 
onwards, LEP increased in energy through the 161 GeV threshold, for the W boson 
pair-production, and reached 209 GeV in 2000. This phase, referred to as LEP2, 
continued the precision measurements, enlarging our understanding of the Standard 
Model, probing the theory to a high level of accuracy. 

The LEP programme has thus set a milestone in the experimental determination 
of the parameters of the SM and has confirmed the theoretical predictions with ex- 
tremely high precision. Once the detectors and the accelerator have been dismantled, 
the LEP tunnel will house the Large Hadron Collider (LHC), a proton-proton accel- 
erator with -y/s = 14 TeV, which is expected to begin operation in 2007. This chapter 
describes briefly the LEP collider and, in more detail, the ALEPH detector and its 
performance. 



2.1 The LEP collider 
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Figure 2.1: An overview of the LEP accelerator. From Ref. [p^]. 



2.1 The LEP collider 

The Large Electron Positron collider (LEP) lay between 100 and 150 m underground on the 
Swiss- French border, near Geneva. It was a 26.7 km long storage ring designed to accelerate 
electrons and positrons and produce collisions at four points where general-purpose detectors 
were placed to observe and record the resulting new particles. These detectors were ALEPH, 



DELPHI, OPAL and L3. Figure 2.1 shows the LEP accelerator complex along with the detec- 



tors. For a detailed description of LEP see Refs. |3^, |3^ , here only an overview is given. 

Beams of counter-rotating electrons and positrons were injected into LEP at an energy 
of 22 GeV from the Super Proton Synchrotron (SPS). Once the particles were in the highly 
evacuated beam-pipe, they were accelerated by superconducting NbTi radio frequency (RF) 
cavities and normal Cu cavities operating at a frequency of 352 MHz. To maintain a circular 
orbit, the beams were bent by 3368 dipole magnets and focused by 816 quadrupole and sex- 
tupole magnets. When particles follow a circular trajectory, they loose some of their energy 
by emitting synchrotron radiation. Since this effect is inversely proportional to the square of 
the accelerator radius, it is necessary to build very large accelerators. Hence the large radius 
of LEP helped to minimise the amount of energy that cavities had to give back to the particles 
to attain the nominal ^/s. 

The normal mode of operation consisted of four equally spaced bunches of electrons and 
positrons which crossed each other every 89 /is (the collision frequency was thus / ~ 11 kHz). 
With approximately 10^^ particles per bunch, typical total beam current values were 6 mA, 
or 750 /iA per bunch. At the interaction points (IPs), where collisions are produced and the 
dimensions of the bunches have to be as small as possible to enhance the collision rate, each 
bunch was approximately 1cm long, 200 /xm wide horizontally {6x) and 8 fim vertically [Sy). 
This yielded a typical LEP2 instantaneous luminosity of 10^^cm~^s~^. The instantaneous 
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luminosity parameter contains all the accelerator capabilites and is proportional to the rate of 
interactions (collisions). It is defined as: 

where A^bunch is the number of bunches and A^c± is the number of electrons/positrons per bunch. 
The luminosity was measured independently by LEP and the four experiments with dedicated 
subdetectors. 



2.2 The ALEPH detector 

The ALEPH (Apparatus for LEp PHysics) detector was designed to cover as much solid angle 
as possible around the point where the beams were made to collide, with high granularity and 
hermeticity. This is to ensure that all particles emerging from the e"'"e~ collisions are measured 
and identified. In practice, a coverage of ~3.9srad was achieved, allowing the detailed study 
of Standard Model physics at LEP and searches for new physics. Its shape was cylindrical, 
built around the interaction point, with 12 m in diameter and 12 m in length, weighing over 
3000 tons and having some 700 000 readout channels. 

Particles originating at the IP and traveling outwards traversed through several subde- 
tectors arranged in cylindrical layers. They first passed through a series of three low-density 
tracking devices, designed to measure the trajectory of charged particles, and then encountered 
high-density calorimeters where all but muons and neutrinos were completely stopped deposit- 
ing all their energy. Only muons could penetrate to the last subdetector, the muon chambers, 
being therefore tagged as such. The tracking subdetectors were immersed in a highly uniform 
1.5 T magnetic field created by a 6.4 m long and 5.3 m diameter superconducting solenoid. The 
trajectories of charged particles are curved in the strong magnetic field, describing a helix which 
spirals around the beam axis. This curvature provides a measurement of their momentum. 



Figure shows all the subdetectors forming the ALEPH detector. Moving radially away 
from the beam-pipe, there is the Microvertex Detector (VDET), the Inner Tracking Chamber 
(ITC) and the Time Projection Chamber forming the tracking volume of the detector. They 
are followed by the Electromagnetic Calorimeter (ECAL) and then the magnet so as not to 
degrade the ECAL energy resolution. The iron return yoke for the magnet was segmented 
into layers and instrumented with streamer tubes acting as a Hadron Calorimeter (HCAL). 
Finally, surrounding the HCAL, two additional double layers of streamer tubes formed the 
muon chambers. Another three pairs of detectors LCAL, SiCAL and BCAL provided low 
angle coverage and measured Bhabha scattered electrons and positrons. They were located 
very close to the beam-pipe, BCAL giving measurements of instantaneous luminosity from 
outside the detector, while LCAL and SiCAL measured integrated luminosity. The ALEPH 



detector is described in Refs. p9| , ^ pj] and its performance is described in Ref. . 

The coordinate system is defined as follows: the origin is the geometrical centre of the 
detector, also the nominal IP; the x-axis is horizontal and points to the centre of the LEP ring; 
the z-axis is along the electron beam direction and makes an angle of 3.59 mrad upwards with 
the local horizontal; the y-axis then makes an angle of 3.59 mrad with the local vertical. This 
small displacement with respect to the local vertical is due to the fact that the LEP tunnel 
was tilted by 1.4% for geological reasons. Given the detector and event geometry, it is often 
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more useful to use spherical {T,6,(f)) or cylindrical coordinates (p,0,z), defined as: 
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(2.2) 



where the polar angle 6 is measured with respect to the +2;-axis and extends between and vr. 
The azimuthal angle is measured in the {x, y) plane starting at the x-axis. The coordinate 



system is depicted in Fig. 2.3 



2.3 Charged Particle Tracking 
2.3.1 Vertex Detector 

The Vertex Detector (VDET) was the innermost subdetector of ALEPH. It was very 
close to the IP to provide highly accurate measurements of charged particles' trajectories. It 
specially served to identify very short-lived particles, such as i?-mesons, made of 6-quarks, 
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which travel a very small distance before decaying (typically ~400 fim). Its principal aim was 
to determine whether a particle came from the IP or rather originated from a displaced vertex. 
It was upgraded in 1995 to prepare for LEP2. With respect to the original VDET, it had 
improved radiation hardness, covered a wider angle and introduced less passive material. The 
main physics purpose of the upgrade was to improve the efficiency of 6-tagging. At LEP2 
energies, the Higgs boson is expected to decay mainly to 6-quarks (BR(if — > bb) > 70%), 
hence 6-tagging is crucial in Higgs boson searches. 

The final VDET consisted of two approximately cylindrical layers. The inner layer, resting 
on the beam-pipe at 6.3 cm from the beam, was made of 9 'faces', each face measuring 40 cm in 
length along the 2;-axis and formed by six double-sided silicon microstrip detectors, or 'wafers'. 
The outer layer, at 11.0 cm from the beam, had 15 faces and ensured a maximal lever-arm 
between the beam-pipe and the Inner Tracking Chamber. Figure p.4| shows a full view of the 
VDET and the geometry of the faces from an end-on view. 
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Figure 2.4: Left, full view of VDET. Right, end view showing the position of the faces Ell 
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The two sides of each wafer had readout strips arranged orthogonally to each other. When 
a charged particle passed through the wafer, a VDET hit was formed giving a two dimensional 
measurement: the r0 side provided a measurement while the z side recorded z. The readout 
pitch on the r0 side was 50 /im and 100 /xm on z. However, the spatial resolution was improved 
by using the signal recorded on several strips to determine the position of each hit down 
to 10 /im in r0 and 15 /xm in z. This precision allowed very precise measurements of particle 
lifetimes, as short as 300 fs. The angular acceptance was such that a particle with | cos^^| < 0.95 
hit at least one layer of the VDET. 



2.3.2 Inner Tracking Chamber 



The ITC was a 2m long cylindrical multiwire drift chamber. There were eight concentric 
layers of wires running parallel to the beam [£) direction. Its inner radius was r = 16 cm, where 
it supported the outer layer of the VDET, and its outer radius extended to r = 26 cm. Its 
full volume was filled with a mixture of argon (80%) and carbon dioxide (20%) at atmospheric 
pressure. Each of the 960 sense wires, arranged in eight layers, was held at a positive potential 
of ~2kV and was surrounded by five field wires and a calibration wire all held at ground 
potential. This formed an hexagonal 'drift cell'. The cell structure is shown in Fig. 2.5. 
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Figure 2.5: Drift cell geometry of the ITC @. 



O Sense Wire 
* Field Wire 
© Calibration Wire 
Calibration 'zigzag' 



When a charged particle traversed an ITC cell it ionised the gas and produced electrons 
which were attracted to the sense wires. The ionisation charge drifted to the sense wire with an 
average velocity across the cell of 50/im/ns. This negative current pulse was detected at both 
ends of the sense wire. The z coordinate was thus obtained by the difference in arrival time 
of the signal at each end. The r0 coordinate was deduced by converting the drift time into 
a drift distance using a parametrisation of the non-linear relationship between the two. The 
spatial resolution was around 150 /im in the r0 coordinate and only 5 cm in the z coordinate. 
Each wire could only produce one hit per event, and therefore there were a maximum of eight 
hits per track. 

The drift cells were small so that they could provide information extremely quickly: the 
cell size determines the maximum drift time. This was achieved at the cost of having many 
wires, readout channels and the consequent problem of fitting more electronics. But then, the 
ITC was able to provide the only tracking information for the first-level trigger (see Sec. |2.5| ) 
in less than 3 /xs after a beam crossing. 
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2.3.3 Time Projection Chamber 

The Time Projection Chamber (TPC) |^ was the main tracking detector in ALEPH providing 
three-dimensional points for track reconstruction. It was a 4.4 m long cylindrical drift chamber 
with its axis parallel to the beam, an inner radius of 31cm and an outer radius of 180 cm. 
There were three main components: the field cage, made up by two cylinders (inner and 
outer), two circular end-plates, and eight 'feet' which supported the TPC weight, attaching it 
to the magnet. A diagram of the TPC is shown in Fig. p.6| . 




The chamber was divided into two halves by a central membrane kept at high potential 
(— 27kV). This membrane, along with the cathode planes of the end-caps (held at a potential 
near ground) and the inner and outer field cages, produced a uniform electric field in the ±z 
direction. The passage of a charged particle ionised the gas (ArCH4) and the electron cloud 
drifted towards the nearest end-cap following the 11 kV/m electric field. The parallel magnetic 
field within the tracking volume ensured that drift electrons did not diffuse radially and had 
tight helix trajectories (the clouds were contained in the r0 plane). 

The TPC end-plates were formed by three consecutive wire grids. See Fig. 2/7 for a diagram 
of an end-plate edge. The first one was the gating grid which will be discussed shortly. Next 
was the cathode grid, which was grounded and was responsible for the drift of the ionisation 
cloud. And finally, electrons were collected on the sense grid. It consisted of alternating sense 
wires, held at ground potential, and field wires, held at a high positive potential, which formed 
conventional proportional multiwire chambers with cathode 'pads' to read out the pulses. 
There were 21 rows of cathode pads precisely positioned behind the wires and capacitively 
coupled to them. Hence the z coordinate was calculated from the arrival time of the signal 
on the pad and the known drift velocity of electrons in the drift volume. The drift velocity 
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Figure 2.7: Schematic view of a TPC sector edge (left) showing wire grids, pad plane and potential 
strips. On the right, the gating grid of the TPC showing electric field lines in the (a) open and (b) 
close gate states. From Ref. |4C]. 



was constantly measured by laser calibration. The r and coordinates were obtained from 
the position of the pads involved in the measurement. Hits could overlap in r0, e.g. if two 
particles ionised the same r0 region at different 2, producing two pulses in the same pad. The 
pad recognised the different pulses by their different time profile allowing separation of the hits. 
Thus the TPC could measure up to 21 three-dimensional points for a track at polar angle 9 of 
90° with a resolution of 180 /xm in rcj) and 800 /im in z. The average resolution decreased with 
particle momentum and 9. Events with low momentum and low polar angle, like two-photon 
processes, have the worst resolution. 

If an ionisation avalanche was produced close to the sense wires, positive ions could drift 
back to the main TPC volume and distort the local electric field leading to track distortion. 
The gating grid in the end-caps was used to prevent charge build-up in the drift volume of the 
TPC. It was thus located between the drift volume and the cathode wires so that if the gating 
grid was in the 'closed' state, positive ions did not enter the drift volume. In the 'open' state, a 
negative potential was applied to the gating wires and the gate was transparent to the passage 
of drifting charged particles. In the closed state, alternate gate wires were kept at positive and 
negative potentials creating dipole fields which prevented the passage of charged particles. See 
Fig. for a diagram of the gate electrostatic configuration in both open and closed states. 
The gate was open 2 /is before a beam-crossing and stayed open for 45 /is if the Level-1 trigger 
reached a 'yes' decision (Sec. ^I5|) . This allowed the ionisation electrons to be collected. A 
negative decision from the Level-1 trigger closed the gate until the next beam-crossing. 

Furthermore, the TPC wires were also used to measure the specific ionisation of charged 
tracks since the pulse height on the sense wires is proportional to the specific ionisation energy 
loss (d-E/dx). The coordinates {r,(f),z) along a particle trajectory provide a measurement of 
its momentum. The dE/dx for each particle depends only on its velocity in a certain material, 
so by combining the dE/dx and the momentum measurements the mass of the particle can 
be inferred and thus its identity derived. There were a total of 338 possible measurements 
of dE/dx for tracks within the angular acceptance of the TPC. A typical resolution of 7.2% 
was achieved in hadronic events. Not all tracks had dE/dx information since a minimum of 
50 wire samples were required to obtain a reasonably accurate measurement. Thus particles 
with polar angles near or vr did not have identification information from the TPC. Figure p.8| 
shows how the dE/dx information from the TPC is used to obtain particle identification. 
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Figure 2.8: On the left, the measured dE/dx as a function of particle momentum for a sample of about 
40 000 tracks. Lines are the fitted parametrisation to a modified Bethe-Bloch formula for electrons, 
muons, pions, kaons and protons. On the right, the average separation in standard deviations between 
different particle types as a function of momentum. It has been computed using all tracks in hadronic 
Z decays which have at least 50 dE/dx measurements. From Ref. |42]. 



2.4 Calorimetry and Luminosity Monitors 
2.4.1 Electromagnetic Calorimeter 

The main purpose of the Electromagnetic Calorimeter (ECAL) was to measure the energy 
of photons, electrons and positrons in an event. The ECAL covered polar angles of | cos^| < 
0.98 and was capable of measuring tracks throughout the polar angle range of the TPC. It was 
inside the superconducting solenoid to avoid a degradation of its performance resulting from 
particles decaying in the uninstrumented coil. 

The ECAL was composed of a 4.7 m long barrel and two end-caps, all with similar prop- 
erties. The calorimeter consisted of 45 layers of lead with a proportional wire chamber in 
between each layer. The signals from the wire chambers were read out by cathode planes ar- 
ranged into towers pointing to the interaction point. There were more than 73 000 such towers, 
each subdivided into three 'storeys' which each yielded an energy measurement. 

An incoming electron, positron or photon interacted in one of the 2 mm lead sheets and 
generated an electromagnetic shower by bremsstrahlung in the case of electrons and by e^e~ 
pair-production in the case of photons. When the particles from the shower traversed the 
wire chambers, they ionised the gas and the cathode pads collected the pulses from the wires. 
This allowed a position measurement of the shower development to be made from the readout 
of the cathode pads. The information from the wires with their faster readout was used to 
make a trigger decision. The choice of lead as absorber was due to its short radiation length 
(Xo = 6.4g/cm2 H). The total thickness of lead was about 22Xo, or ~40 cm, which ensured 
that over 98% of the energy of a 50 GeV electron was contained in the lead sheets. Hadrons 
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and muons traversed the ECAL and could deposit some energy inside it since it was about one 
interaction length thick. 

The large number of towers provided high granularity which was needed for spatial separa- 
tion of the electromagnetic showers and helped to improve particle identification. Since each 
tower gave three energy measurements, the shower profile could also be studied, adding more 
information to discriminate between electromagnetic and hadronic showers. An ECAL layer is 



illustrated in Fig. 2.9. 
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Figure 2.9: Detail of an ECAL layer. There are 45 such layers Q. 
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The end-caps met the barrel at polar angles of ~40°. There were a number of cables in this 
region which reduced the total depth in radiation lengths. The energy resolution was degraded 
by around 30% in this region as a result. 



2.4.2 Hadron Calorimeter 



The Hadronic Calorimeter (HCAL) was used to finally stop and measure the energy of 
hadrons as well as tag the trajectories of muons. Its structure was very similar to that of the 
ECAL, comprised of a barrel and two end-caps. However, it was slightly rotated with respect 
to the ECAL to avoid an overlap between the inactive 'cracks' of both detectors. This can be 



seen in Fig. |2.10| , which offers a view of both ECAL and HCAL. 

The HCAL was made of 23 iron layers separated by plastic streamer tubes. It weighed in 
excess of 2500 tons and particles at normal incidence encountered 1.2 m of iron, equivalent to 
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Figure 2.10: The ALEPH calorimeters. The HCAL end-caps measure about 9 m in diameter |4C]. 



7.2 interaction lengths for hadrons. The iron acted as the return yoke of the magnetic flux 
from the solenoid and also as the main mechanical support for the ALEPH detector. 

A hadron reaching the iron layers interacted through the strong force with the nucleons 
of iron and decayed into a shower of low energy particles. Hadronic showers are broader and 
deeper than electromagnetic showers which explains why the HCAL needed to be thicker than 
the ECAL. The ionisation caused by the shower in the tubes was collected by cathode pads 
which were grouped into 4068 projective towers pointing to the IP, like the ECAL. The pad 
electrodes gave a measurement of the integrated energy from each tower, but did not provide 
any information on the r coordinate, the depth of the shower, which is very important to 
determine the hadronic shower profile. A digitised cathode strip running along each streamer 
tube provided the necessary r0 measurements by indicating which plane of the HCAL had 
recorded a hit. This enabled a better discrimination between muons and hadrons by their 
different hit pattern (see Sec. |2.6.3| ). 



For pions at normal incidence the energy resolution is parametrised by |^ 

^ ^ 0.85 
E GeV) 
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2.4.3 Muon Chambers 

Muons traversed both calorimeters with a characteristic trail of hits with no shower develop- 
ment. Surrounding the HCAL there were two further double- layers of streamer tubes known 
as muon chambers. Only muons could penetrate that far after the HCAL and these chambers 
acted as a muon tracking detector. The two layers were separated by more than 40 cm to 
enable a measurement of the track angle with an accuracy of ~15mrad. 

2.4.4 Luminosity Monitors 

The luminosity delivered by LEP was determined by using three pairs of specific detectors. An 
accurate measurement of the luminosity L is necessary in order to obtain precise reaction cross 
sections from the event count rates. Cross sections, giving the probability of an interaction to 
occur, are usually given in barns, defined as 1 b = 10~^'^m~^. 

These detectors measured the number of small angle Bhabha events, coming from the QED 
process e+e~ — > 7* — >■ e+e", whose cross section is known very accurately, to calculate the 
luminosity^. 

For N events and cross section a: 



where s in GeV is the centre-of-mass energy and 6'min/max are the angular bounds which define 
the acceptance cone subtended by the monitors. 

Each luminosity monitor consisted of two units on either side of the interaction point to 
detect the outgoing back-to-back e~ and e"^. LCAL and SiCAL were inside ALEPH itself, both 
at about 2.5 m away from the IP. BCAL, with 12 layers of tungsten and scintillator counters, 
was about 7 m away, outside ALEPH. 

SiCAL had 12 annular tungsten-sheets sandwiched between silicon detectors. The active 
region covers polar angles from 24 to 58mrad. LCAL was similar in construction to the 
ECAL — lead sheets interspersed with multiwire proportional chambers — , covering angles 
just beyond SiCAL, from 45 to 190 mrad. LCAL can be thought of as a low angle extension of 
the ECAL. LCAL and SiCAL had low event rates so they provided an integrated luminosity 
measurement. However, for LEP2, additional background shielding was required close to the 
beam-pipe which reduced the angular acceptance of SiCAL to 34 mrad. So, at high energies, 
SiCAL lacked enough data to provide an accurate measurement of / Cdt and was used only 
to extend the calorimeter hermeticity in ALEPH. LCAL was therefore the only source of 
integrated luminosity. 

BCAL, on the contrary, having a much lower angular acceptance (6'min ~ 7 mrad) due to 
its position, could detect a larger number of Bhabha events and had high enough statistical 
precision to provide an instantaneous luminosity value. BCAL was therefore used as an online 
luminosity monitor of the background conditions. 



^The same process, via exchange of a Z boson, introduces some corrections in the total cross section which 
were poorly known at the start of LEP. By restricting the measurements to low 6, the contribution from Z 
production is effectively reduced, and the rate of events is much larger. This is why these detectors sit very 
close to the beam-pipe. 




with 




(2.6) 
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2.5 Trigger System 

At LEP, there was approximately one beam crossing every 89 fis but the detectors could not 
possibly record all those events at such an enormous rate. A decision on whether the event 
was 'interesting' was required to activate the readout system and record and store the event. 
However, of all the beam crossings, most resulted in 'uninteresting' events, like beam inter- 
actions with gas molecules in the accelerator, or off-momentum beam particles scattered by 
a collimator close to the IP. Only a small fraction of the beam crossings actually consisted 
of electron-positron interactions. The trigger system was designed to accept these events and 
reject the others. It also reduced the dead-time of the detector that results from readout. 

The trigger system was built in three stages, the first two were hardwired logic to achieve 
a very fast response and the third one implemented in off-line software. Each one based their 
'yes' or 'no' decision on different subdetectors depending on their capability to produce fast 



information about an event. Table |2.1| summarises the subdetectors used and the decision 
time. 



Stage 


Decision Time 


Max. Rate (Hz) 


Information Used 


LEP beam crossing 


~ 11000 




Level- 1 


5 /is 


~ 10 


Hit Patterns in ITC 
Pad/ wire readout 
from ECAL+HCAL+LCAL 


Level-2 


50 fis 


~ 5 


Hit Patterns in TPC 


Level-3 


62 ms 


1-3 


All subdetectors 



Table 2.1: Summary of the ALEPH trigger system. 



The first level trigger (Level-1) delivered a positive decision within 5 /is, much less than the 
beams crossing, if a good charged track had been registered by the ITC and/or the calorimeters 
had an energy deposit. This positive decision allowed the TPC to collect the drift electrons 
and partially reconstruct the tracks. If the Level-1 gave a 'no' decision, the TPC gate grid was 
closed and all the subdetectors reset for the next beam crossing, avoiding any dead time. The 
Level-2 needed around 50 fis to form a decision based on the TPC tracks: it required good 
tracks (exceeding a hit threshold) pointing to the interaction region. If the event satisfied the 
Level-2 conditions, the whole detector output was read out. 

The Level-1 and Level-2 criteria valid for the 1999 and 2000 run periods is the following, 
with only minor modifications on the threshold values for previous years: 

• Single Electron/Photon Triggers. Drift chamber (ITC for Level-1, TPC for Level-2) track 
segment pointing at an energy deposit in the ECAL of 200 MeV for electrons or a 2 GeV 
energy deposit for photons (with no track). 

• Total Energy Trigger. ECAL energy deposits should exceed certain thresholds: 5.5 GeV 
in the barrel or 4.5 GeV in either end-cap or 1.7 GeV in both end-caps. 

• Back-to-Back Track Trigger. Two ITC (or TPC for Level-2) tracks pointing in opposite 
directions from the IP. 



• Single Muon Trigger. ITC or TPC track with at least 4 HCAL hits. 
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• Bhabha Luminosity Triggers. Hits in the LCAL with opposite azimuth. 

In fact, the trigger achieved an efficiency greater than 99% in accepting the 'interesting' events, 
such as Z — i> gg, Z — > , Bhabha scattering, or two-photon events. 

2.6 Event reconstruction 

Event reconstruction is the process by which raw data from the different subdetectors is trans- 
formed and correlated to form physically meaningful final data. It mainly consists of track 
reconstruction, which assigns each hit in the tracking chambers to a particle track, and energy 
flow analysis, which associates charged particle tracks with energy deposits in the calorimeters, 
enabling particle identification and improving the overall energy resolution. Event reconstruc- 
tion begins just after the data is stored on tape and is run by JULIA (Job to Understand Lep 
Interactions in Aleph) |^ . Particle identification, as well as the rest of the user analysis of the 
events is performed by the computer program ALPHA (ALeph PHysics Analysis package) [^] . 

2.6.1 Track reconstruction 

The information from all three tracking detectors is used to reconstruct a charged particle 
trajectory, or track. The reconstruction algorithm starts in the TPC. It links nearby hits to 
form segments. Several segments will be connected to form a track if the hypothesis that a 
helix is described is satisfied. This TPC track will then be extrapolated down to the ITC and 
VDET where matching hits will be assigned to the track. This solves any ambiguities that 
could arise from the two-dimensional (r0 and z) hits in the VDET. 

A final track fit is then performed taking into account the segment errors by means of 
a Kalman filter technique |^ to provide a fully reconstructed track. Table |2.2| shows the 
momentum resolution achieved by incorporating the diferent subdetector segments into the 
track for a Z — > fj^~^fi~ sample of events. 



Tracking Detectors 


Max. Hits 


a,Jp^, (GeV/c)-i 


TPC 
TPC+ITC 
TPC+ITC-fVDET 


21 
21+8 
21+8+2 


1.2x10-3 
0.8x10-3 
0.6x10-3 



Table 2.2: Transverse momentum resoultion cTp^/pl for pt in GeV/c when the information from the 
different subdetectors is incorporated into the track reconstruction p^]. The number of maximum 
hits available for each subdetector corresponds to a track with 90° polar angle. 

A 'good' track, as applied to the analysis described here, is defined as follows: 

• A minimum of four hits in the TPC, which restricts good tracks to the range in polar 
angle: | cos 6*1 < 0.96 

• Transverse impact parameter: |(io| < 1cm 

• Longitudinal impact parameter: |2;o| < 5 cm 

• Transverse momentum: pt > 2MeV/c 
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2.6.2 Energy Flow 

The energy flow algorithms are designed to provide an accurate calculation of the total visible 
energy and momentum of an event. The energy resolution is fundamental in SUSY searches 
since their characteristic signature is large missing energy or momentum. 

The event is first of all 'cleaned' by removing possible cosmics (the muons from cosmic 
showers initiated in the upper atmosphere) and electronic noise. Clusters of energy deposition 
are built in the calorimeters by adding together 'storeys' (in the ECAL) or towers (in the 
HCAL) which have energies above 30 MeV and share a common edge or corner. Good tracks, 
as defined in the previous section, are then extrapolated to these calorimeter clusters. This 
forms the so-called energy flow objects, which are groups of topologically connected tracks and 
clusters. Several such objects are defined: 

• Electrons. They are identified using the ECAL shower profile. Their energy is calculated 
using their track momenta and mass. If the associated calorimeter cluster energy is 
greater than 3 times the expected energy, that extra energy is associated with a photon 
from bremsstrahlung. 

• Muons. Identified by 400 MeV energy deposits per HCAL plane and up to 1 GeV in an 
ECAL cluster. 

• Photons. Neutral energy clusters in the ECAL. If the invariant mass of two photons 
matches the mass of neutral pion is assumed. 

• Charged and Neutral Hadrons. The remaining charged tracks are attributed to charged 
hadrons, and the remaining energy to neutral hadrons. 

Finally, this analysis provides a list of objects which is expected to be an accurate repre- 
sentation of the true particles involved. The kinematic variables calculated from energy flow 
objects provide the basis for any further analysis of the event. For example, the visible mass of 
an event (Mvis) is the invariant mass of all energy flow objects in that event; or the transverse 
momentum (pt) is the transverse component of the vector sum of momenta of all energy flow 
objects. 

The energy flow resolution as a function of the polar angle is shown by the lower set of 



markers in Fig. |2.11| , calculated for hadronic events Z qq. The distribution of the recon- 
structed energy is fitted to a Gaussian giving an average relative resolution of ~7%, reaching 
a maximum of 15% for cos 6' > 0.98. The ALEPH hermeticity is also shown in that figure, as a 
measurement of the reconstructed energy normalised to ^i. The gradual decrease with cos 6* 
in the measured energy for hadronic events arises from the high multiplicity of such events. 
Events with cos 9 > 0.9 are more likely to have objects going down the beam-pipe, leading 
to a decrease in the measured energy. Otherwise, Bhabha and muon-pair events have low 
multiplicity and are very well reconstructed. However, there is a sudden drop in the measured 
energy for Bhabha events at cos 6* ~ 0.98. This corresponds to the insensitive region where the 
LCAL meets the ECAL end-cap. 

The high hermeticity of the detector and the precise measurements of the energy contained 
in the events, permit searches for new physics to be carried out successfully in ALEPH. 



2.6.3 Lepton Identification 



Lepton identification has a major role in the analysis described here, and needs to be described 
in some detail. 
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Figure 2.11: Energy flow resolution normalised to one as a function of the polar angle 9 {\ = 71/2 — 9). 
An average energy flow resolution of 7% is achieved. The three sets of markers at the top of the 
plot show the reconstructed energy as a fraction of the total input energy {^/s). The sharp drop at 
cos 9 ~ 0.98 for Bhabha events is due to the gap region between the LCAL and ECAL. From Ref. [p9|. 



Electron Identification 

Electron identification, as described in the previous section, relies on the energy profile de- 
posited in the ECAL and additionally on the di?/dx information from the TPC Three quan- 
tities, called estimators, are defined such that the measured value of a variable (g) is compared 
to the mean expected value ((q')), taking the resolution (cXg) into account. The estimator (Q) 
will have a standard Gaussian distribution: Q = centered at zero and with width unity. 

The estimators are: the transverse shower shape Rt-, based on the compactness of the shower; 
the longitudinal shower shape Rl, and i?/, based on the ionisation energy loss in the TPC. 

The electron shower in the ECAL tends to be more compact in the transverse direction 
than that of a hadron like the pion. This leads to the definition of the Rt estimator, based 
on the compactness of the shower. The momentum p of the electron track is known from the 
TPC data. That track is extrapolated to the ECAL and the energy in the four towers closest 
to the extrapolated point is summed, giving E^^. The energy fraction E^^/p is the variable used 
in the Rt estimator: 

^ Ei/p{me&s) - Ei/p{exp) 

The second estimator -Rl is defined to take into account the longitudinal development of 
the electron shower, using the inverse of the mean depth of the energy deposition in the ECAL 

^ ^ Ea ^ ^ XL(meas) - Xijexp) 

' Eti EU 
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where the energy is summed over i storeys and j segments, S is the mean depth of energy 
deposition on that segment. 

The estimator Rj is similarly defined, by comparing the measured dE/dx to what is ex- 
pected for a particle with the mass and charge of an electron. 

It is also important not to misidentify electrons produced originally in the IP with photon 
conversions, i.e. photons which interact with the VDET, ITC or TPC material and produce 
an e"'"e~ pair. This e"'"e~ pair does not have to be symmetric, and sometimes only one particle 
is energetic enough to be registered. Since there are many SM backgrounds involving photons, 
it is important not to take a conversion for a 'genuine' electron coming from the signal SUSY 
process. This is achieved by trying to identify pairs of tracks coming from a single displaced 
vertex. To tag conversions it is required that: 

• the distance between the two tracks in xy plane at the point where they are parallel and 
closest together, is less than 2 cm, 

• the distance between the two tracks in the z direction at that same point is less than 
3 cm, 

• the invariant mass of the two tracks at that point assuming they are both electrons is 
less than 40 MeV. 

Finally, for an object to be identified as an electron it must pass the following criteria: 

• It must be a 'good' track and not compatible with a photon conversion 

• p>2 GeV/c 

• i?T > -3 

• \Rl\ < 3 

• \Ri\ < 3 

The average efficiency of the procedure has been evaluated to be 65% in hadronic events . 
See Figs. pj.l2| a and b for the dependence of the efficiency with the transverse momentum and 
the polar angle. No evidence for a strong correlation between the two is found. 



Muon Identification 

Muons have a different hit pattern in the HCAL than hadrons; they only interact with HCAL 
material by ionisation, not nuclear interactions, following a very narrow path, rather than 
extensive transverse areas. Muons are also expected to reach the outer regions of the HCAL, 
which means that planes further away from the IP are more likely to be fired by a muon than 
a hadron. Finally, the track extrapolated from the TPC, the narrow hits in the HCAL and 
the possible hits in the muon chambers aid in the correct identification of tracks as muons. 
Furthermore, these patterns are independent of the muon momentum for muons in the range 
5 — 50 GeV/c. All muons above 3 GeV/c are able to penetrate out to the muon chambers. 

Tracks are extrapolated from the TPC assuming they are muons defining a 'road' through 
the HCAL with a width of three times the uncertainty due to multiple scattering. A plane 
in the HCAL is expected to fire if the road intersects it. A plane is said to have fired if a 
digital hit lies within the width of the road. A hit in the muon chambers is assigned to the 
extrapolated track if it lies within four times the uncertainty due to multiple scattering. The 
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Figure 2.12: The efficiencies of the electron identification procedure as a function of (a) the transverse 
momentum and (b) the polar angle. The electron sample used comes from real data at the Z peak 
selecting photon conversions with high purity, (c) The efficiency of muon identification using both the 
HCAL and muon chambers information as function of the polar angle. The sample used is again real 
data Z IJ-^fJ-^ where high purity is obtained by selecting muons with TPC and ECAL information 
only. From Ref. ||]. 



depth of penetration of the track can be assessed by looking at the last ten planes, out of the 
total expected to be fired, that actually register a hit. Thus, a track can be identified as a 
muon if it satisfies the following conditions: 

• iVfi,e/iVe,p > 0.4 

• A^e^p > 10 

• A^io > 4 

• -^mult < 1-5 

• At least one associated hit in the muon chambers 

where A^fire, A'exp, A^'io and Xmuit are, respectively, the number of actually firing planes, the 
number of expected planes, the number of firing planes within the last ten expected and the 
average hit multiplicity per fired plane. 

Applying these constraints on a dimuon events sample, an average efficiency of 86% is 
found. Figure |2.12| c shows the muon identification efficiency as a function of the polar angle. 



The two dips are due to the magnet supporting legs, which are not instrumented with muon 
chambers. 

Tau Identification 

The tau particle, the heaviest known lepton, can decay both leptonically or hadronically via 
the exchange of a W boson. This happens with a mean lifetime of 290 fs, or 87 /im, just enough 



to be measured at ALEPH for boosted taus. The main decay modes are listed in Tab. |2.3| . If 
the final state consists of one charged particle it is referred to as one-prong. This is the most 
likely decay mode with 85% probability and gives some relatively clean signal in the leptonic 
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Mode 


X 


Main Decay 


BR(%) 






T~ —>■ e^Z/cZ^r 


18 








17 


V^Ilc-piUUg, 


h 


71^ Ur 


13 


h + n' 




25 






T" 7l~27l^Ur 


10 




/i + 37r° 




1.5 




?>h 




9.5 


3-prong 


3/i + (> IvrO) 


7r^7r"'"7r~7r°z/T- 


5 


5-prong 


hh + (> 07r°) 


T- ?>h-2h+v^ 


0.1 



Table 2.3: Tau lepton main decay modes with the corresponding branching ratio (BR), rounded to 
nearest integer. X is the final state for r — > Xvr- h stands for charged hadron and in general means 
a TT^, although is also possible. Neutral pions decay promptly to two photons and charged pions 
are detected in the HCAL. The one-prong decay modes add up to 85%, the three-prong to 15% and 
the five-prong to ^^0.1%. Adapted from Ref. (5|. 



mode. On the other hand, for Lorentz-boosted taus decaying in the three-prong mode, it is 
usually very difficult to reconstruct the final particles (tt^ and tt" in general) due to the small 
opening angle between the three. 

Tau identification is performed by clustering the energy flow particles into two or four jets, 
depending on the analysis, with the Durham algorithm [Q. A jet will be a tau-jet candidate 
if it contains one or three good tracks (or two if it contains an identified electron, to allow for 
asymmetric photon conversions) and if the jet invariant mass is less than 2 GeV. A cut on the 
cluster-radius as a function of the momentum is also available. 

It is important to define well the structure of jets, be it in Higgs, W or r decays. Several 
algorithms exist that operate by joining adjacent energy flow objects to form a jet. The Durham 
jet clustering algorithm used here, described in Ref. [Q, uses the invariant mass Mij of two 



objects i and j as a proximity measure. The procedure is as follows: 

1. all possible pairs between particles are formed and the quantity 

M2 = 2[mm{Ei, Ej)]\l - cos %) (2.9) 

is evaluated, where Ei and Ej are the energies of the objects and 6ij is the angle between 
them. 

2. the pair with the minimum value of M^j is replaced by a pseudo-particle with energy 
Ei + Ej and momentum pi + pj 

3. steps 1 and 2 are iterated over the set of pseudo-particles and remaining objects until 

Vij = > Vent (2.10) 

vis 



Thus the parameter ?/cut sets the proximity condition for two tracks to be grouped together. 
The higher the value of ?/cut, the easier it is to form jets with the available particles: so a larger 
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7/cut will give fewer jets. The definition of Mij is such that while the separation 6ij is important, 
the association of two objects is still possible if the value min(ii^j, Ej) is sufficiently low. 

For instance, if we want to force the event into four jets, the value ?/cut for which the event 
is grouped from three into four jets receives the name of j/34. A genuine four jet event will 
generally have a large value of 7/34, thus tightening the value of 1/34 will reject processes that 
do not generally produce four jets. 



2.7 Event Simulation 

Any analysis in searches for new particles requires two event samples to be generated. One is 
the signal process we are interested in, for example slepton pair-production, and the other is the 
set of all known Standard Model processes that could mimic the signal process in the detector, 
the background. Finally, of course, there is the actual data collected with the detector which 
will decide if the signal has been found in the data and is distinguishable from background. 
Monte Carlo (MC) event samples for signal and background are generated with specific physics 
generators and processed in ALEPH in four stages. 

1. The physics generator emulates e'*"e~ collisions and produces the final state particles and 
four-momenta of the process or processes in question. The output is compliant with 
a set of kinematic rules imposed by KINGAL, an ALEPH program that will generate 
the kinematics of all particles in the final state to simplify further processing. This 
information is then passed onto the next stage. 



GALEPH, the ALEPH detector simulation program, based on GEANT3 which 
simulates the response of the ALEPH detector to the final state particles. GALEPH 
contains a detailed and updated description of the geometry and composition of the 
detector and simulates the interactions of the generated particles with the sub-detectors 
on their passage through the different materials. The output format of GALEPH is the 
same as a real event. 



As for real data, the next stage is the reconstruction program, JULIA which fa- 
cilitates the later analysis by forming physical objects. Unlike data, the reconstructed 
Monte Carlo files contain information on the true simulated event, thus enabling useful 
comparisons to be performed. 

The final stage is the analysis itself where, based on the process under study, a set 
of variables is calculated for each event to try and make it distinguishable from the 
background and provide any necessary measurement. This is done with ALPHA, the 
ALeph PHysics Analysis package [H^. 



Usually, background files are generated and made available for the collaboration. In the 
analysis described here, background files were not generated by the author. 



Chapter 3 

Searches for six-lepton final state 
topologies in GMSB 

The fi NLSP scenario is favoured in GMSB models with the number of messen- 
ger famihes greater than one. In this case, pairs of selectrons and smuons may be 
produced at LEP2 and decay immediately to final states with two leptons and four 
taus and missing energy. This type of decay depends crucially on the mass differ- 
ence between the sleptons and staus, but no less important is the neutralino mass, 
which enters the process as a mediating particle either on shell (two-body decays: 
^ ^ ^TTi) or off-shell (three-body decays: Zr — Zrfi). Both processes are 
studied in this chapter. Up to now, only evidence for two-body decays has been 
looked for at colliders. The analysis described in this thesis is the first experimental 
search for the three-body process. 

This final state topology is interesting because it covers a very sensitive area of 
parameter space, specially in the low gravitino mass range, i.e. for negligible stau 
lifetimes. Other searches, for cascade decays of neutralinos or direct slepton decays, 
are unable to cover this region of parameter space and thus the appeal to see if Nature 
has realised SUSY in this region. 

After an introduction to the decays and the motivation for the search, the experimen- 
tal searches are described in detail and since no evidence for this process has been 
found in the highest energy ALEPH data, limits on the cross section and sparticle 
masses are set. 
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Figure 3.1: The decay of selections and smuons when the hghtest stau is the NLSP, can occur via a 
two-body decay to an on-sheh neutrahno, or via a three-body decay directly to the NLSP through a 
virtual neutralino. Both cases are addressed in this work. 



3 . 1 Phenomenology 

In the context of Gauge Mediated SUSY Breaking (GMSB) models, final states with six leptons 
may appear when the lightest stau fi is the NLSP and pairs of lightest selectrons e^ or smuons 



/iR are produced. The pair-production of sleptons at LEP is shown in Fig. [L^a. The decay 
chain under study in this thesis is shown in Fig. |3.1| . 

The final state is thus two leptons I, either electrons or muons with opposite sign, four taus, 
and missing energy carried away by the undetected gravitinos. In the case of a light gravitino, 
rriQ < lOeV/c^, the NLSP will decay promptly after production. The analysis presented here 
will be restricted to this case. If the gravitino mass is between 10 and around 500eV/c^, the 
fi will present some sizeable decay length in the detector and the topology would then have 
kinks or large impact parameters from the fi decays. This analysis 'with lifetime' has been 
performed by OPAL in Ref. . 



The decay of the selectron or smuon /r, occurs via a neutralino. It has usually been assumed 
in experimental searches that the neutralino would be on-shell with a mass range between: 



m 



I — nil > ITT'x > ^fi + 



where, / stands for e or fi. And hence the decay would follow two-body kinematics: 
Two — body : e^e~ —>■ IrIr ^ Ix ^ ~^ IttQ IttG 



(3.1) 



(3.2) 



But specially important is the fact that this decay can proceed through an off-shell neu- 
tralino with the special kinematics of three-body decays: 



> mi — mi > mf^ + mr 



(3.3) 



Three — body : e"'"e — > /r ^ Itti Itti IttG IttG (3.4) 
Three-body decays of selectrons or smuons into the lepton, tau and stau have been theo- 



retically studied in Ref. 54 1 and this thesis is the first experimental search for this particular 
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Figure 3.2: Pair-production cross sections at ^/s = 207 GeV for the lightest neutralino, selectron 
and smuon in the (a) two-body and (b) three-body scenarios. Selectrons (dark shaded area) may 
have larger production cross sections than neutralinos (light shaded area) in the two-body case, but 
never in the three-body case. The lightest stau production cross section, not relevant for six-lepton 
topologies, is very similar to the smuon cross section since they do not have t-channel contributions. 



type of decay. In this case, there is the advantage of eliminating the upper constraint on the 
neutralino mass (left hand side in Eq. which can now lie well above the /r, mass or the 
beam energy at LEP. This search is therefore sensitive to an unexplored range of neutralino 
masses at LEP2. 



Figure ^]2| shows the production cross section as a function of particle mass for neutralinos, 
selectrons and smuons in the two-body and three-body scenarios. The /ij^/ij^ production cross 
section is only larger than that of e^e^ for masses greater than ~90GeV/c^ in the three- 
body scenario. So one may expect more fi~^ fi~ tt f if i events than e"'"e~rrfifi with three-body 
kinematics, although this is not guaranteed. Smuons are only produced via s-channel exchange 
of a 7 or Z, while selectron and neutralino pair-productions are model dependent via the 
additional t-channel exchange of a neutralino and a selectron respectively, see Fig. |1.6| . The t- 
channel contribution to the neutralino cross section is constructive, hence the generally higher 
neutralino cross section in the three-body case with respect to the two-body case, where the 
exchanged selectron is heavier than the neutralino. On the contrary, the interference between 
s and t diagrams is destructive in the selectron production, so that light neutralinos cause the 
t-channel to dominate and reduce the total cross section. 

As a result, three-body /r decays are disfavoured at the production level with respect to 
neutralino production with a slepton NLSP. The decay of selectrons and smuons with a stau 
NLSP has to compete in the available parameter space with direct slepton NLSP production 
{i —>■ iG) and neutralino production with a slepton NLSP (x ^ HG). However, the six- 
lepton final state topology offers a complementary channel for the discovery of supersymmetry 
in GMSB models and is sensitive to exclusive areas of parameter space. A discussion on the 
available parameter space in GMSB models for this type of topology follows. 
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3.1.1 Available parameter space 

In GMSB models with N5 > 1 the fi NLSP scenario dominates. This can be clearly seen in 
Fig. p.3| , where Region 1, corresponding to a x NLSP is only dominant for N5 < 2. Regions 2 
and 4 correspond to a fi NLSP scenario, where the six-lepton final state topology is possible. 
Three-body decays of /r have more parameter space available (Region 4) than two-body decays 
(Region 2) for N5 > 1. It is also noticeable how Region 4 extends down to relatively low values 
of tan /3 ~ 5 independently of A. For even lower values of tan (3, the /r — fi mass difference 
becomes less than the tau mass and the three sleptons coexist as co-NLSP (Region 3). 

In the stau NLSP scenario, if the mass ordering is mj^ —mi> m^, then the two-body decay 
~^ dominates over /r — > IG. The latter decay is suppressed by the gravitational nature of 
the gravitino. If the mass order is > —mi> rrif^ + rrir the three-body modes /r —>■ Wfi 
will dominate^ as shown in Fig. p.4| . The decay to /r — ^ Ix is impossible in this region of 
masses, and only /r IG may become available in pathological cases where the slepton mass 
is very close to the stau mass and the decay to the stau is thus not possible. Only then the 
direct decays from sleptons to gravitinos may be realised overcoming the very low gravitational 
coupling. Once tan f3 is greater than ~10 the stau mass becomes much lighter than the slepton 
masses due to mixing effects and the three-body decay modes are restored. It is found that 
all slepton decays proceed through /r — > Wfi and Zr IG is avoided if the difference between 
the slepton and the stau masses is greater than 5GeV/c^. 



3.1.2 Kinematics of two- and three-body decays 

In the two-body scenario (Eq. the mediating neutralino decays independently into t^t^ 
or T^f^ with equal probability. One expects therefore no asymmetry in the charges of the 
two staus: there will be equal numbers of like-signed staus and opposite-signed staus. But in 
the three-body scenario, it is important to distinguish between the different charge channels, 
since the decay widths r(Zj^ — > Z~r+ff ) and — > l~T~f^) can actually be quite different. 
This asymmetry is caused by the different coupling of the neutralino to the stau depending 
on the stau mixing. If the lightest stau is purely right-handed, there is no asymmetry and 
the largest hypercharge states and are predominantly produced in equal amounts as 
primary taus. But if the left-handed stau content of the lightest stau increases, i.e. the mixing 
angle Of increases, then the wino (W3) nature of the neutralino will prefer the state with more 
TjJ content and the asymmetry will be created. 

Assuming a fully bino neutralino, which is generally almost true in GMSB models, and a 
near threshold production (small /r momenta), one can deduce the following relations^: 



r(/R - i-r-f^) 



(m2 -m? )2 



■ n ^X n 

sm Uf H — COS tif 

mr 



(3.5) 



Til, 



■r) 



m 



2 \2 



cos Of + 



m^ 



sin Of 



m, 

'R 



(3.6) 



Both the 'slepton-charge flipping' channel (Eq. |3.5| ) and the 'slepton-charge preserving' 
decay (Eq. |3]^) are suppressed for increasing neutralino mass, with other parameters held 



^There could also be three-body decays Fr — > viPtTi through ofF-shell charginos xt^ although they are 
strongly suppressed due to the large chargino mass and the very small coupling of to I'lX^- 
^After some algebra from the calculations in Rcf. p3] 
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Figure 3.3: In the plane (A,tan/3) and for N5 = 1,2,3,4 (a, b, c and d respectively), regions of 
X NLSP (1), i co-NLSP (3) and n NLSP with < m^-^ (2) and m 



X > (4)- 



Mm has been 



fixed at 3A and /i is positive. The heavy solid lines denote the boundaries between these regions. 
Regions 2 and 4 are those relevant for the six-lepton final state topology, where it has to compete 
with X ~^ 'T^i s-iid direct production processes. The cross-hatched region at high tan (3 is excluded 
because electroweak symmetry is not correctly broken. The shaded regions are excluded because 
rrif < 76GeV/c^ (vertical shading), < 84GeV/c^ (horizontal shading) or rrih < 95GeV/c^ (or 
niA < 83GeV/c^) (diagonal shading). The dot-dashed contours are where the chargino mass is 95, 
200 or 350GeV/c^, while the dotted line is the contour of rue^ = 90GeV/c^. From Ref. |^]. 
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Figure 3.4: Slepton branching fractions as a function of tan /? in the three-body scenario for two 



ranges in Am 
(or mv 



■ mr 



Only for highly virtual neutralinos (Am^^- > 30 GeV/c^) and tan /? < 10 



< 5GeV/c ) may the channel /r IG be open and dominate. Otherwise the three- 



body decay modes dominate. The asymmetry between I r f^*") and r(/j 



l-T+f; 



is 



explained in Sec. 3.1.2. The condition > mj^ — mi > mf^ + mr is only possible for tan /? > 5 and 
Ns > 2. 



fixed. Figure shows the asymmetry between the two decay modes as a function of tan/3, 
which is roughly inversely proportional to the stau mixing angle Of, and as a function of 
Am^j- = — mj^. It is clear from the plots that: 



(3.7) 



so that same-sign t'^t'^t^t^ states are suppressed compared to opposite-sign states t^t^t^t^ . 
This effect will be most significant for large m^/mj^ and for small negative values of cos^f. 
For mj^ — TJif^ < lOGeV/c^, | cos^f| usually ranges between 0.1 and 0.3 (small mixings). 

The kinematics of the event will ultimately depend upon the mass differences: mj^ — irtf^ , 
\'fni^ — m^l and — m^^, neglecting lepton masses. Figure ^]5| shows the energy distributions 
for the lepton and the primary tau in both two- and three-body decays with m^-^ — nif^ fixed 
at ~llGeV/c^. In the two-body scenario, as one would expect, the energy of the primary 
lepton is fully dependent on the /r — x mass difference whereas the energy of the tau lepton 
depends on — mf^. But in the three-body scenario, these energy distributions prove to be 
completely independent of the neutralino mass: i.e., no matter how virtual (m^ — Tni^) the 
neutralino is, the lepton and primary tau energies are determined only by the total available 
energy -^^beam - mf. 

The analysis described here will assume negligible stau lifetime (ttIq < lOeV/c^). Could 
the sleptons nevertheless have a long lifetime? The probability that a slepton Zr with energy 
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Figure 3.5: Energy distributions, with arbitrary vertical axis, for lepton (electron or muon) and 
primary tau at the generation level. The dotted histograms represent the two-body decay modes 
~^ a-iid the full histograms, the three-body decays /r — > Wfi. Signal points are denoted by: 
mi^/m^/mf^. In both cases Am = m^-^ — ruf^ is fixed at ~llGeV/c^ and the neutralino mass is 
changed to see the effect. 



E wiU travel a distance x in the lab frame before decaying is given by: 

/I rneV\ / F2 

P[x) = e-^/^ with A = CT-fp = 0.2 — - 1 mm (3.8) 




It can be shown that the decay length for sleptons produced at LEP2 (where E is the 
beam energy, and thus the Lorentz factor is not much greater than 1) is only measurable if 
Am = mj^ — TUf^ is less than one GeV. It is therefore assumed in this analysis, without loss of 
generality, that the sleptons will not have a sizeable decay length. 

Examples of how selectron and smuon production with six-lepton final states might look 
are shown in Fig. |3T 



3.2 Data sample and Monte Carlo simulation 

Six-lepton events with missing energy in GMSB scenarios constitute the events of interest in 
this analysis, and they are called 'signal' events. Once they have been described and the new 
physics behind them is understood, it is necessary to know how would they look were they 
to be produced in a collision in the ALEPH detector. This is performed by means of Monte 
Carlo simulation. The same technique is utilised in generating Standard Model processes that 
could mimic or fake the signal events we are interested in. These processes, arising from known 
physics, are called 'background' events. A 'selection' procedure to distinguish between the two 
will be devised such that it retains the maximum number of signal events and rejects as many 
unwanted background events as possible. Finally, the selection is applied to the real data 
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(a) Two selectrons were produced with 99GeV/c^ which then decayed to two electrons (tracks 2 and 5) and two 
neutralinos with mass 90GeV/c-^. Of the four produced taus, three underwent a one-prong decay into pions (tracks 
1, 3 and 9) and one decayed into three pions (tracks 4, 7 and 8). 



DALI_F2 ECM=208 Pch=83.3 Efl=86.3 Ewi=29.6 Eha=26.2 dtrec_08 Run=1999 Evt=712 

ALEPH Nch=10 EV1 = EV2 = EV3 = ThT=0 Detb- 71 




-300 cm Z 300 cm 

Y'=cos (210) *Y-sin(210)*X 



(b) Two smuons decaying via a three-body process. The two corresponding muons are tracks 1 and 2. One of the 
primary taus decayed into a tt"*" (track 3) and the other into three pions (tracks 5, 6 and 7). The secondary taus 
coming from the staus decayed into a muon (track 4) and into three pions (tracks 8, 9 and 10). There is a photon 
EncfcotlverSBte in the inner wall of the TPC, between tracks 3 and 4. 

Figure 3.6: Signal Monte Carlo events at ^/s = 208 GeV. The detector is shown in the 'fish-eye' view 
to magnify the tracking relative to the calorimeters. 
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collected by ALEPH. A measurement of the degree of evidence for the signal process under 
investigation can then be obtained by comparing the selected number of events in the data 
and the expected number of events from the background simulation. 

3.2.1 Data sample 

The data analysed in this thesis were collected with the ALEPH detector between 1998 and 
2000 with a total integrated luminosity of 628 pb~^ at centre-of-mass energies between 189 
and 209 GeV. Table. p.l| gives the detailed luminosity at each energy. The LEP accelerator 
was pushed to its limits during year 2000 and ALEPH collected 11.6pb~^ of data between 200 
and 204 GeV, 197.3 pb"^ between 205 and 207 GeV and 7.3 pb"^ at more than 208 GeV. The 
excellent performance of both the accelerator and the detector must be acknowledged. 



Year 


1998 


1999 


2000 


(v/i) (GeV) 


188.6 


191.6 


195.5 


199.5 


201.6 


205.0 


206.7 


JCdt (pb-i) 


173.6 


28.9 


79.9 


87.0 


44.4 


79.5 


134.3 



Table 3.1: Average centre-of-mass energy and corresponding luminosities of the analysed data sample 
for data collected by the ALEPH detector from 1998 to 2000. 



3.2.2 Signal 

The ALEPH SUSY Task Force generally uses SUSYGEN pBI as the signal Monte Carlo pro- 



gram for the generation of supersymmetric events in the simulation. All searches described 
in this thesis were performed using SUSYGEN, except the six-lepton search described in this 
chapter. SUSYGEN is unable to produce the type of events where selectrons and smuons are 



heavier than the stau and the decay is a three-body mode. Instead, ISAJET 7.51 fS^] which has 
the full capability to generate /r — * /rfi events was used. ISAJET had to be interfaced with 
KINGAL to allow the generated events to be processed by the standard ALEPH reconstruction 
routines. 

In the GMSB framework, the minimum number of parameters to describe a complete 
sparticle spectrum is six: A, Mm, (or v^Fq), N5, tan/3 and sign(/i). But phenomenologically 
there are only three relevant parameters in the search for slepton decays: the masses of the 
SUSY particles involved, itlj^, and m^^, which will completely determine the kinematics 
of the event. Instead of an artificial grid with values for the three masses, the scan over the 
six theoretical parameters described in Sec. ^]l| was used to obtain the interesting values for 
the masses. Thus having scanned millions of possible GMSB models, the chosen signal points 
were obtained by imposing the appropriate hierarchy in masses (Eqs. p.l| and |3.3| ) and the zero 
lifetime condition < lOeV/c^ for the stau NLSP. 

Some reduction was needed in the number of selected signal points, and this was performed 
in the three mass-planes involved (m^, J, (m^, mj^) and (m^-^^, J by requiring a separation 
between the points of more than 2GeV/c^ in any one plane. A high density of points was 
achieved without unnecessary overlapping. Only stau NLSP masses greater than 70GeV/c^ 
were used. In Fig. the signal points are ordered with decreasing m^^ in the horizontal axis 
and the values of their three defining masses are shown in the vertical axis to illustrate the fine 
coverage of all possible kinematical regions. 
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Signal Point Signal Point 



Figure 3.7: The signal points in GMSB parameter space were obtained from the scan described in 
Sec. 5.1. The relevant kinematical parameters, the mass of sleptons, neutralino and stau, are shown 
here for each signal point in the two-body scenario (left) and three-body scenario (right). In the three 
body scenario Am = m^^ — nif-^ > 5 GeV/c^ to ensure that no slepton decays into IG. 



A total of 385 and 457 points in the minimal GMSB parameter space were used as signal for 
the two-body and the three-body scenarios respectively. For each point, a total of 1000 events 
were produced at y/s = 189, 206 and 208 GeV. Taking into account the beam energy limit in 
the production of sleptons, this represents a total of 2 million reconstructed signal events. 

Finally, the simulated signal points in the two- and three- body scenarios are shown in 
Fig. in the (m^,mfj) plane. Also shown are the present limits derived from other GMSB 
searches (described in Chapter ^). This plot demonstrates the reach of three-body slepton 
decays over an area in parameter space where the neutralino is heavier than the beam energy 
and is not excluded by existing limits. This is the main motivation to use this topology and 
cover that region of parameter space. In the case of two-body decays, one is limited by the 
kinematical conditions and is further constrained by the suppressed slepton production cross 
section with respect to neutralino production. 



3.2.3 Background 

Standard Model processes have to be simulated and reconstructed in ALEPH in order to 
compare data and 'known physics' predictions; only when the first shows disagreement with 
the latter does the possibility for new physics arise. It is therefore very important to understand 
and be able to simulate correctly all possible processes that could resemble the one sought. 



Table p.2| lists for each centre-of-mass energy the cross section, number of reconstructed events 
and MC generator for each of the analysed backgrounds. The background generation was 
performed by other members of the ALEPH collaboration, corresponding to at least 20 times 
the collected luminosity in the data. 

Cross sections for some of these processes at LEP2 energies are shown in Fig. p.9| . For 
comparison, typical SUSY particle cross sections are ~0.1 pb or less if the particles are produced 
at threshold (m ~ y/s/2), as can be seen in Fig. |3l2 . 



50 



Searches for six-lepton final state topologies in GMSB 




Figure 3.8: Signal events were generated at each one of the two-body (left) and three-body (right) 
points displayed here in the neutralino-stau mass plane. The hatched areas mark the present exclusion 



from GMSB searches in that plane for < lOeV/c , as explained in Fig. 5.3 



Two-fermion processes 



These are e"'"e~ ^ Z/7 ^ ff processes, where f is a quark or lepton (Fig. |3.1CI| a). They were the 
dominant background at LEPl when the centre-of-mass energy was the Z mass, and thus the 
cross section was maximal. As LEP energies have increased away from the Z peak and above 
the WW threshold, the dominant background has become four fermion processes. At LEP2, 
e"'"e~ —>■ Z will often be produced with an ISR (initial state radiation) photon at low angle. 
This is the so called radiative return to the Z and leads to missing momentum signatures. 
Only T~^T~ and qq events could resemble the six-lepton final state topology and constitute a 
relevant background to the signal. 

Four-fermion processes 



W pair-production e+e W+W such as shown in Figs. p.lO| c and d, constitutes a source 



for jets and leptons since W can decay into qq (70%) and iu (30%). Similarly, e+e~ — > ZZ* 
(with Z* = Z or 7) is a major background for SUSY processes. The possible decay modes of 
the Z boson qq (70%), (10%) and uu (20%) form final states which overlap with different 
SUSY signals. Also possible are e"^e~ — > Wez/ and e"'"e~ — > Zee, see Figs. |3.10| e and f, which 
both involve photon radiation from one of the initial beam particles. The photon propagator 
contributes to the cross section with weight 

^ (3.9) 



Ebeam-Ee(l " COS 6* 



tag; 



where and ^tag are the scattered electron energy and polar angle. Thus events with low 
values of ^tag will have larger cross sections. The scattered electron is then more likely to 
go undetected down the beam-pipe leaving its characteristic missing transverse momentum. 
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Figure 3.9: Cross sections for Standard Model processes as a function of the centre-of-mass energy. 
Taken from Boudjema and Mele in Ref. |36|. 



Conversely, the neutrino produced in Wez/ processes arises from a virtual W whose propagator 
is given by 

(3.10) 



cos 9,, 



E\yea,m Ey ( 1 

This expression resonates for large neutrino scattering angles. So that characteristic Wev events 
will have large missing momentum and a visible mass close to the W mass. Although the cross 
section for this type of process is around one pb, they constitute a difficult background to 
reject because six-lepton final states can be quite similar. For the same reasons, Zee events are 
typically identified by a very energetic electron, visible mass close to the Z mass and missing 
momentum. 



Two-photon processes 

LEP is not only an e"'"e~ collider but also a high energy photon collider. Beam particles are 
constantly radiating photons which may interact with one another with high probability. The 
cross section for these e+e~ e"'"e~ ff processes is of the order of a nb. Figure |3.10^ shows 
the relevant Feynman diagram. As discussed above, the majority of these radiative events are 
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(a) Two-fermion with ISR 



(b) ZZ* dominant contribution 




+ e 



AAAAAAAAAW- 



AAAAAAAAAW^ 



(c) WW s-channel 



(d) WW i-channel 
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(e) Wez/ 



(f) Zee 
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(g) Two-photon 



Figure 3.10: Feynman diagrams of SM background processes: (a) two-fermion, (b-f) four-fermion 
and (g) two-photon . 
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Process 


189 GeV 


196 GeV 


202 GeV 


206 GeV 


MC 




N 


a(pb) 


N 


a(pb) 


N 


a(pb) 


N 


a(pb) 


Generator 


ee 


2000 


966 


1700 


894 


400 


842 


480 


811 


BHWIDE 




75 


8.27 


50 


7.57 


50 


7.08 


50 


6.80 


KORALZ 


TT 


75 


8.21 


50 


7.54 


50 


7.04 


50 


6.74 


KORALZ 


qq 


400 


99.40 


600 


89.74 


125 


84.54 


125 


80.60 


KORALZ 


WW 


108 


16.59 


100 


7.61 


100 


17.46 


100 


17.54 


KORALW 


Wei/ 


15 


0.66 


20 


0.75 


20 


0.83 


20 


0.88 


PYTHIA 


zz 


40 


2.76 


50 


2.85 


50 


2.84 


50 


2.81 


PYTHIA 


Zee 


450 


99.11 


250 


99.04 


210 


98.94 


400 


98.72 


PYTHIA 


77 — > ee 


1200 


3800 


2000 


3840 


2000 


3970 


600 


190"f 


PHOT02 


77 nfj, 


1293 


3550 


2000 


3620 


2000 


3660 


600 


185t 


PHOT02 


77 TT 


210 


431 


300 


441 


500 


454 


300 


90.00t 


PHOT02 


77 qq 






520 


585 


606 


592 


602 


630 


PHOT02 



Table 3.2: Standard Model background processes analysed in the six-lepton topologies search. At 
each centre-of-mass energy, the number of produced events N in units of 10^ and the corresponding 
cross section is given for each process. In the 'tagged' two-photon event generation it is required that 
the visible mass is greater than 2.5GeV/c^ and the scattering polar angle greater than 5mrad. The 
77 — > i£ subsample at ^/s = 206 GeV was required to have visible masses greater than lOGeV/c^, 
hence the lower cross sections (marked with t). The MC generators are BHWIDE |5|], KORALZ pj, 
KORALW H, PYTHIA pi and PHOT02 [||]. 



produced when the two final state e"'"e~ continue through the beam-pipe and are not detected 
(referred to as 'untagged' events). They can be distinguished by their low visible mass and 
particle multiplicity. Nevertheless if one of the final state electrons does appear in the detector, 
generally at very low angles, these events can fake SUSY signals with large missing energy. 



3.3 Signal Selection 

Two analyses have been designed to search for /r — * Ix and /r — > /rfi in the fi NLSP scenario. 
No independent search has been made for either two- or three-body decay modes since both 
decays ultimately depend on the mass difference between the produced sleptons and the stau: 
Am = mr — rrif, . 

The case of large Am, when Am > lOGeV/c^, is the most favourable since the primary 
leptons are energetic, multiplicity is somewhat larger and these events contain large visible 
mass. When Am is small. Am < 10 GeV/c^, there are fewer identified leptons, greater missing 
momentum and much less visible mass. 



3.3.1 Description of the selections 

The final state topology consists of four taus, two of them very energetic, and two leptons 
(electrons or muons) which might be soft if the mass difference between the sleptons and the 
stau is small. Only in the large Am case is it possible to reject events containing fewer than 
two identified leptons, which is a very successful cut against two-fermion backgrounds. In the 
case of small Am a more detailed study of the event is required and a procedure was devised 
in which the event is clustered into four tau jets using the Durham algorithm. 
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Six-lepton Selection for large Am 



Charged tracks 
Preselection Identified leptons 
Visible mass 



4 < Ar,h < 12 
Ni>2 

M^is < O.SSy/i 



Anti-77 



Visible mass 
Missing pt 



Mvis > 25 
£/t > 0.02v/i 



Anti-4f 



Hadronic mass 

Neutral hadronic energy 



M,i,„2^ < 120 

Ehado < o.i7i5;ch 



Anti-2f,ZZ Jet clustering 



Anti-WW ({qiv kinematics 



7/23 > 2 X 10-3 
^34 > 1 X 10-^ 



Xww > 4.5 



Table 3.3: Selection criteria for six-lepton final states with large Am 
GeV/c^, momenta in GeV/c. 



m 



In 



Masses are in 



The selection criteria for the large Am case are summarised in Tab. 3.3 and for small Am in 



Tab. |3]J. The principal rejected background is cited beside the cut, although there are usually 
various backgrounds affected simultaneously. 

The variables used in both the small and large Am analyses are described below. 

• The number of charged tracks A^ch and good tracks A'good 

The number of charged tracks, also known as event multiplicity, peaks at around 20 for 
four-fermion events. Signal events have much lower multiplicity. The number of good 
tracks (as defined in Sec. 2.6. 1| ) restricts N^h to tracks with a minimum of four TPC hits. 



small impact parameter and pt > 2MeV/c. 
Number of identified leptons Ni 

Total number of electrons and muons in the event. Lepton identification was described 



in Sec. |2.6.3| . This is a very sensitive parameter to Am. In the large Am case it helps in 
rejecting Wez/, ZZ, 77, and two-fermion processes which all peak at zero Ni if taus are 
involved. 

Variables M^is, p{, E^h, ^hadO> ^had, and Mhad-2£ 

The visible mass M^is is calculated by adding the energies and momenta of all energy fiow 
particles in the event. The missing momentum is identical in magnitude and opposite 
in direction to the total momentum. Its projection onto the transverse plane defines pi- 
Two-photon events have very small values of Mvis and and these variables are very 
effective against such processes. E^h is the contribution to the visible energy coming from 
charged particles, i^had from all hadrons and ^'hadO from neutral hadrons. To cut against 
four-fermion events, the Mhad-2£ variable is defined as the total hadronic mass with the 
energy of the two most energetic leptons subtracted. The normalised Mvis against Mhad 



distribution for the WW background is displayed in Fig. |3.11| a. The small Am signal is 
shown to have much less hadronic energy than four-fermion events and significantly less 
visible mass. 



3.3 Signal Selection 



55 



Six-lepton Selection for small Am 




Good tracks 


2 < A^good < 10 


Visible mass 


M,is < O.lOv^ 




Visible mass 
Missing pt 


Mvis > 0.05^ 
j/t > O.lOv/i 




Hadronic mass 


Mhad < 50 


Anti-4f 


Hadronic mass 


Mhad-2£ < 0.25v^ 




Charged tracks in two tau jets 


Ar-+- < 6 




r jets with good tracks 


jygood > 2 


Anti- TT 


Charged tracks not in tau jets 


A',^h°^^ > 2 




Thrust 


Thrust < 0.95 



Table 3.4: Selection criteria for six-lepton final states with small Am = m^^ — nif^. Energies are in 
GeV/c^, momenta in GeV/c. 



Cut values of the jet clustering algorithm ?/23 and 2/34 

The Durham jet clustering algorithm |50| was described in Sec. 2.6.3| . Here, the sensitive 



variables are 1/23 and 1/34, which are the ?/cut values (from Eq. |2.10| ) necessary to force the 
event from two to three jets and from three to four jets respectively. An example of the 
distribution of the rr and ZZ background along with the signal for the large Am case is 
shown in Fig. |3.11| b in the plane (?/235l/34)- 



WW qq£z/ chi-squared Xww 

The kinematics of these events can be exploited to suppress the background by defining: 
,2 _/'"^qq- "^w Y f mu-mw \^ / p; + -Emiss - 94GeV/c ^ ^ 



" V lOCivT?; \T0GeVj?) ^[ lOGeV/c J ^^'^^^ 

where mqq is the hadronic mass, i.e. the mass of the event after removing the leading 
lepton, rriiiy is the mass of the leading lepton and the missing momentum, pi is the 
momentum of the leading lepton and -Emiss = a/s ~ -^vis is the missing mass. WW events 
are likely to occur at small XwW' ^^"^ therefore be rejected by requiring a minimum 



Xww events to be selected. See Fig. p. 12 



Four-tau jets variables N^^^^ A^s°°<^ and N^f^ 

These variables are defined by forcing the event into four tau jets. A^^"*"^^ is the number 
of charged tracks inside the two most energetic tau jets. It is used against Wez/ and ZZ 
backgrounds. Two very effective anti-rr variables are: A|°°^, which is the number of tau 
jets that contain good tracks; and A"^^"^, the number of charged tracks not associated 
with any of the two most energetic tau jets. 

Thrust 

The thrust is defined as: 

Thrust = max ^^^'^^^ (3.12) 

" Ei \Pi\ 

where the sum runs over all reconstructed particles. The thrust ranges from 0.5 for 'open' 
(spherical) events to one for events with two well coUimated jets. It is helpful in rejecting 
TT events with thrust values very close to one. 
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Figure 3.11: (a) Visible and hadronic mass distributions for WW background (open circles) and 
six-lepton signal points (three-body decay mode m^-^ :m^:mfi = 97:123:94 GeV/c^) in the small Am 
selection after a subset of cuts has been applied, (b) Background distributions of rr (open cir- 
cles) and ZZ (triangles) events in the (2/23, 2/34) plane before the cuts on those variables in the 
large Am selection. The six-lepton signal is shown as empty squares, in a two-body decay mode: 
mj -.m^-.mf^ = 95:87:73 GeV/c^. 



The distributions of total background, data and signal after the preselection and anti-77 
cuts for the large and small Am selections are shown in Figs. p.l2| and p.l3| for ^/s = 206 GeV 
with an integrated luminosity of £ = 126.5 pb~^. 

3.3.2 Optimisation 

As can be seen in Fig. p.l2| , to decide where to place the cut, one has to reach a compromise 
between the number of accepted background events, which should be small, and the efficiency 
in selecting signal events, which should be kept high. Such a compromise can be obtained 
analytically by means of the 'ergs prescription', described in Ref. |6^. 

The number N of expected events from any given process with production cross section a 
if the experiment records C integrated luminosity and the selection efficiency is e, is given by 

= eaC (3.13) 

Nevertheless, in a search experiment there are two independent contributions to the number 
of observed events: the expected background b and the expected signal s. The number of 
'candidate' events n selected in the data after the selection has been applied can be described 



by Poisson statistics |64 



P^n;s + b)= ^^+/^^" e-(-+^) (3.14) 



which gives the probability to observe n events in the data if we expect s + b. If 6 is small 



3.3 Signal Selection 



57 




Figure 3.12: After applying the preselection and anti-77 cuts, the distributions of the remaining 
variables used in the large Am selection. The data (dots) at ^/s = 206 GeV are compared to the 
background Monte Carlo (filled histograms). The dashed histograms show the signal distributions 
(for the values m^-^ : : m^j = 95 : 87 : 73GeV/c^) in an arbitrary normalisation. The location 
of the cut is indicated with an arrow. The cut values were optimised and fixed before analysing the 
data, as described in Sec. [3.3.2 . 
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Figure 3.13: After applying the preselection and anti-77 cuts, the distributions of the remaining 
variables used in the small Am selection. The data (dots) at ^/s = 206 GeV are compared to the 
background Monte Carlo (filled histograms). The dashed histograms show the signal distributions 
(for the values mj^ : : m^j = 97 : 123 : 94GeV/c^) in an arbitrary normalisation. The location 
of the cut is indicated with an arrow. The cut values were optimised and fixed before analysing the 
data, as described in Sec. [3.3.2 . 
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and n is compatible with this expectation, then it is reasonable to assume that the observed 
candidate events come from background and not from signal. Thus an upper limit on the 
number of possible signal events can then be calculated. Usually the upper limit a„ for n 
candidates observed, is given with a 95% confidence level (C.L.), by solving: 



k=0 



"'ra 



0.05 



(3.15) 



This is the probability to see n or less candidates if we expect a„. The coefficients a„ = 
3.00, 4.74, 6.30 . . . are the solutions to Eq. |3.15| with n = 0, 1, 2, . . . For example, if no candidate 
events are observed in the data, every signal leading to 3.00 or more expected events is excluded 
at 95% CL. This follows from the fact that the number of observed events from a signal process 
described by a Poisson distribution with a mean of three, has a 5% chance of fluctuating down 
to zero. 

The selection procedure should be optimised to yield the signal efficiency and the back- 
ground expectation corresponding to the best obtainable limit. As regards the background 
expectation, the position x of the cut on the discriminating variable being optimised is deter- 
mined by minimising Ng^. This is the expectation value of a„ in the absence of any signal, i.e. 
if all candidates are background: 



-bix) 



n=0 



b{x) 



-b{x) 



3.00 + 4.746(x) + 6.30 



h\x) 



+ 7.75 



h^{x) 



+ 



(3.16) 



where x is the position of the cut to be optimised, h{x) is the expected number of background 
events surviving the cut and n is the number of observed events. The number h{x) is determined 
after all other cuts on the selection have been applied to the total background Monte Carlo 
sample. Equation. p.l6| is just the sum of the coefficients a„, weighted by the Poisson probability 



to observe n events if h background events are expected. 

However, if the number of expected background is large it is useful to recalculate the 
coefficients a„ taking into account the background. This procedure is referred to as 'subtracting 
the background', and generally only applies to well modelled and understood backgrounds. It 
is possible to only subtract one type of background events and not the total, i.e. only WW 
expected events or even only 80% of them if one is cautious about the systematic uncertainties 
in the estimation. In this case, if hguh is the number of 'subtracted' background events, Eq |3.15| 
is modified to become 0: 



P„(a„ + h, 



sub) 



0.05 



.... (3.17) 

n \'-'sub ) 

The expected limits when background is subtracted are stronger than the limits produced 
when background is not subtracted. Fig. p.l4| shows the different iVgs as a function of the 
expected background events, if background is fully subtracted or not. The background level 
in the selections developed in this chapter is low and there is no need to perform background 
subtraction. 

Finally, the best position x for the cut on the variable under optimisation is such that the 
expected excluded cross section: 



0"95 



[x)C 



(3.18) 



is minimal. While the numerator depends only on the number of accepted background events. 
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Figure 3.14: The expected limit on the signal events as a function of the expected background b, if 
the background is subtracted, and if no background subtraction is performed. 



the denominator only depends (cumulatively) on the number of selected signal events. Thus the 
best possible limit is calculated at the Monte Carlo level without using the data. Figure |3.15| 
shows this process on two selected variables. First, the cut values are set approximately 
by looking at background and signal distributions as shown in Figs. |3.12| and |3.13| . Once the 
chosen cuts give a good background rejection, the list of cuts is frozen and each one is optimised 
sequentially applying the others. 




Figure 3.15: The optimisation procedure for the neutral hadronic energy in the small Am selection 
(left) and the thrust in the large Am selection (right). The signal distribution is shown on the upper 
third, the background in the middle and the minimised CJ95 at the bottom. All other cuts have been 
applied and the arrow indicates the optimal position of the cut on these variables. 



3.3 Signal Selection 
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3.3.3 Selection efficiencies 

The main systematic uncertainties on the efficiencies come from the number of generated events 
in the simulated samples (~3%). The total systematic uncertainty introduced by lepton identi- 
fication variables has been estimated to be < 2% Selected data events and background MC 
were compared to see any difference in the performance of the lepton identification algorithms. 
The estimators were shifted by the difference between data and MC, and new efficiencies were 
calculated to see the impact. These effects have been taken into account by conservatively 
reducing the selection efficiency by one standard deviation. 

The generated signal data contained selectron and smuon events with the appropriate cross 
sections and branching ratios calculated by ISAJET. Thus no separation between selectron 
and smuon events in the selection efficiency is possible. The efficiency values are then for the 
topology as a function of Am rather than independent measurements of the flavour of the 
event. 

Signal events were only generated at 189, 206 and 208 GeV, thus an interpolation for 
intermediate energies was needed. A check on some signal efficiencies was performed at these 
intermediate energies to parametrise the efficiency as a function of ^/s.lt was found that a 
linear interpolation between centre-of-mass values was appropriate since the calculated values 
were well within the systematic errors on the efficiencies of the generated points. Furthermore, 
only a few points had variations larger than 10% between the efficiency at 189 and 208 GeV. 
Finally, the selection efficiencies are presented in Fig. |3.16| as a function of Am = m^^ — 



Tl 



for both two- and three-body signal points. Efficiencies greater than 30% are assured over the 
full range in Am. The spread in the bands is due to the range in slepton masses generated from 
the kinematical limit down to 70GeV/c^. The spread in the large Am selection when applied 
over two-body signal points arises from small or large slepton-neutralino mass differences. 
It is noticeable here that the same selection applied over three-body signal points shows no 
dependence on the virtuality of the neutralino. 
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Figure 3.16: The efficiency of the small and large Am selections as a function of Am for all analysed 
signal points. 
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6£ small Am 


6i lar^ 


re Am 




/ Cdt 


data 


bkg 


data 


bkg 


189 


173.6 





0.73 





0.25 


192 


28.9 


1 


0.12 


1 


0.04 


196 


79.9 





0.34 





0.12 


200 


87.0 


1 


0.37 





0.13 


202 


44.4 


1 


0.14 





0.23 


205 


79.5 





0.34 





0.12 


206 


126.5 





0.54 





0.18 


208 


7.8 





0.03 





0.01 


TOT: 


627.6 


3 


2.60 


1 


1.07 



Table 3.5: After all cuts have been applied to the data and the background Monte Carlo the numbers 
of events for each selection at the different centre-of-mass energies. The P-values are: P{k > 3;b = 
2.60, s = 0) = 0.48 and P{k >l;b = 1.07, s = 0) = 0.65. 



3.4 Results 

In this section the above selections are applied to the data and the background MC samples. 
The selected candidate events are briefly discussed and the significance of the observed results 
is described. Limits on the production cross section for selectrons and smuons are derived and 
conservative limits on the masses set. 



3.4.1 Events selected in the data 

The number of observed candidates in the data along with the SM background expectation 



at each centre-of-mass energy is presented in Tab. p.5| . Generated Monte Carlo samples for 
background processes were only available at the energies detailed in Tab. p.2| . To obtain an 
expected background value for intermediate energies, the two closest available values were 
scaled to the luminosity of the desired energy and linearly interpolated between them. 

In the 628 pb~^ of ALEPH data between 189 and 209 GeV, the search for six-lepton topolo- 
gies with small Am selects three events when 2.6 are expected and the search for large Am 
finds one candidate event with 1.1 expected. What can one then say about the signal process 
under study? To quantify how likely it is to have signal production in the data, i.e. s 7^ in 
Eq. |3.14| , one can calculate the probability to find n events or more if background only was 
expected: 



00 ,^ n—l ,^ 

P{k >n;b,s = 0) = Y: ^f' = 1 - E fc! (3-19) 

k=n k=0 

This is sometimes referred to as the confidence level (CLf,) or P-value. 1-CLf, measures the 
compatibility of the observation with the background hypothesis. As a convention, it is agreed 
that one has discovered a signal if 1-CL;, is less than ~6 x 10^''. That is to say, the probability 
to observe such an unexpected result if the only source is the 'known' background, would have 
to be ~10~^ to be convinced that the signal is present in the data. The median expectation 
for pure background is l-CLb=0.5. Therefore if CLf, is smaller than 0.5 there is a deficit in the 
background expectation. The above equation yields a P-value of 0.48 and 0.65 for the small 
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Run 


Event 


SM interpretation 


191.6 


50210 


5122 


Zee 


199.3 


51940 


8969 


WW 


201.6 


52100 


3592 


Zee 


191.6 


49630 


10846 


77qq 



Table 3.6: Candidate events selected by the small Am selection (first three) and the large Am 
selection (bottom). 



and large Am selections, respectively. Thus the data is in perfect agreement with the expected 
value from the Standard ModelQ. 

The composition of the expected background for the small Am selection is as follows: WW 
and 77rr dominate with 40% each, the rest has equal contribution from rr, Weu and Zee. In 
the search for large Am topologies, the dominant background is ZZ with 35%, followed by qq, 
rr and Zee. All other background sources are completely eliminated. 

The four selected data events are all compatible with SM sources as listed in Tab. |3.6| . They 
are displayed in Fig. |3.17] . The event selected by the large Am selection at 192 GeV, shown 
in Fig. |3.17| a, has two identified electrons. One is isolated with momentum 6.5 GeV/c and the 
other has 66GeV/c momentum, with 27 GeV of hadronic energy in a cone of 30° around its 
direction. There is large missing momentum in the z direction, suggesting that one particle 
went through the beam-pipe. It could therefore be a 77 event with two low multiplicity jets, 
one of them with a semileptonic decay. Of the three candidate events selected by the small 
Am selection, the two recorded at 192 and 202 GeV, in Figs. |3.17] b and d, present a very 
energetic electron of 15 and 66 GeV respectively, and large missing longitudinal momentum. 
Both invariant masses are compatible with mz, at 59 and 98GeV/c^ respectively. Therefore 
they are possibly Zee events where one of the electrons goes undetected. Finally, the event 
at 200 GeV shown in Fig. |3.17| c is probably a WW event with a nuclear interaction very close 
to the IP. There is an identified electron with 50 GeV energy and a three-prong tau candidate 
with an invariant mass of 1.9GeV/c^, suggesting WW ez/rz/. The remaining two tracks 
have together an invariant mass of 0.1 GeV/c^. 



3.4.2 Limits on slepton production 

In the absence of any evidence for a signal, it is possible to constrain the parameters of the 
theory that would produce such a signal. However, the statistical nature of the counting 
experiment makes it impossible to rule out with 100% certainty any hypothetical process. Thus 
all limits quoted here will be expressed at 95% confidence level, which should be understood 
as: 'if the experiment were to be performed again, the probability to obtain a result which is 
in as bad or in worse agreement between the expectation and the observation is less than 5%'. 



^ Of course this argument relies heavily on the background estimation, so it should be used quoting estimation 
errors. In the case of a 10% fluctuation in the background expectation, the above results for the P-value would 
still be compatible with the background only hypothesis. 
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Figure 3.17: Candidate events. 



3.4 Results 
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Mass limits 



The ratio 



CL, 



CL 



s+b 



llP{k<n;s + b) n E (s + 6) Vfc! 

i i k=0 



]\P{k < n; b) 



n E e-''6V^! 



(3.20) 



i k=0 



is used to derive a lower limit on the hypothetical slepton mass. CLf, was defined in Eq. |3.19| 
and CLs_|_{, follows from the same equation when the number of expected signal events s is also 
included. CLg^h is thus a measure of the compatibility with the signal+background hypothesis. 
The product over i channels makes reference to the different energy bins where data was taken. 
Thus n, s and b change with the centre-of-mass energy i (see Tab. p.5| ). To perform the 
calculation of the observed CL^ and the average expected without looking at the data, the 
FORTRAN package described in Ref. |^ was used. 

The limits on the masses are set when the confidence limit (CL^) reaches 0.05, i.e. the 
probability that a signal-like experiment is less compatible with the signal than what is observed 
is below 5%. The confidence level for the different selections as a function of the selectron and 



smuon masses is shown in Fig. |3.18 



— observed 




expected 
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Figure 3.18: Lower mass limits for selections and smuons from the two six-lepton selections. The 
large Am selection confidence level CL^ is plotted in the upper row for selectron and smuon masses. 
The small Am selection in the lower plots. Solid lines give the observed limit and dashed lines the 
average expected limit for background. 
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The limits on the slepton masses derived from the above confidence level are the following: 



Large Am: mg > 99 GeV/c^ > 98 GeV/c^ 

Small Am: mg > 94 GeV/c^ m^ > 95 GeV/c^ 



Cross section limits 

The calculation of a model independent cross section limit is somewhat more problematic 
given the span of centre-of-mass energies of the recorded data. To provide a combined limit, 
the respective luminosities must take into account the appropriate phase space for slepton 
production at each energy. The best approach for the combination is to scale the luminosity 
of the lower energy data by the ratio of the cross section at that energy to the cross section at 
^/s = 208 GeV. Thus the upper limit on the production cross section at the 95% C.L. can be 
written as: 

N95 



(^95 = (3.21) 

(t{208) 

where Ng^ = 3.00, 4.74, 6.30 . . . from Eq. |3.15| depends only on the number of observed candi- 
dates when no background subtraction is necessary as is the case here. This expression takes 
into account the variation of the efficiency with ^/s as well as the change in luminosity if the 
considered slepton mass cannot be produced at that centre-of-mass energy. Using the ratio 
of cross sections to perform this adjustment in the phase space, is more accurate than using 
the simple f3^/ s factor expected for sfermion production. Cross sections are calculated taking 
into account radiative corrections and other model dependent factors whereas the f3^/s scale 
is only valid at the Born level and for s-channel production. Clearly, (J95 has a dependence 
on the sparticle mass through the efficiency function and the cross sections. Thus for every 
possible slepton mass a value of Ugs is calculated. If the production cross section for a slepton 
at that mass is greater than the upper limit (J95, that test mass and more generally the model 
that produced that slepton are excluded. Or inversely, one can say that only those models of 
the theory which produce sleptons with a < ag^ are possible and still realisable in Nature. 

Figure |3.19| shows the excluded cross sections as shaded areas determined by the upper limit 
(T95 for the six-lepton selections as a function of the hypothetical slepton mass. The minimal 
possible production cross section is also shown. 

The interpretation of these results in the wider context of GMSB models is left for Chap- 
ter ^, where these limits and the limits imposed by other searches in the GMSB framework are 
analysed in the most general scenario. 
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Excluded at 95% CL 




Figure 3.19: 95% C.L. excluded cross section areas as a function of tlie slepton mass for the small 
Am (light shaded area) and large Am (dark shaded area) selections. The minimum production cross 
section for pairs of selectrons and smuons covering all possible six-lepton topologies is also shown for 
comparison. 



Chapter 4 

GMSB searches and present limits 



Gauge mediated Supersymmetry breaking models could be realised at LEP in the 
form of many different topologies. These arise from just two factors: whether the 
NLSP is the lightest neutralino or the lightest stau, and whether the NLSP lifetime is 
negligible (prompt decays), short (decays inside the detector) or long (decays outside 
the detector). Searches to cover all possible NLSP types and lifetimes, considering 
also indirect production, have been developed by other members of the ALEPH col- 
laboration and are described here. Their results will be used in the next chapter 
to constrain possible GMSB models. The possible final states under study are: two 
energetic photons, non-pointing single photons, two acoplanar leptons, large impact 
parameter leptons, detached slepton decay vertices, heavy stable charged sleptons 
and multi- leptons plus missing energy final states. 

Negligible lifetime Short lifetime Long lifetime 




4.1 Description of GMSB topologies 
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NLSP 


e+e — > 


Decay mode 


Signal 


X = cT'~f(3 


Expected topology 


X 


XX 


X^lG 


7(7)^ 


A < ^dct 
A ~ ^dot 
A > ^dct 


Acoplanar photons 
Non-pointing photon 
[Indirect search] 


E 






A < idet 
A ~ idet 


Leptons and photons 
Leptons and displaced 7 


or fi 


U 






A < ^dct 

A ~ ^dct 
A > idct 


Acoplanar leptons 
Kinks and large do 
Heavy charged particles 


XX 


X^d^ UG 


££'{££') ^ 


A < £det 
A ~ ^det 


Four leptons 
Kinks and large do 




IrIr 


/r /rfi 


tt{Utt) ^ 


A < £dct 


Six leptons 



Table 4.1: All possible final state topologies according to the NLSP scenario and decay length A. 
Particles in brackets may be soft or even undetected. Only ^^77 signals were not studied in the 
189 - 209 GeV ALEPH data set. 



4.1 Description of GMSB topologies 

There is a breadth of topologies to study at LEP2 to cover all possible GMSB scenarios if one 



considers not only the NLSP type but also its lifetime |3^. Table |4d] lists all possible searches 
in the different NLSP scenarios.^ 

Depending on the gravitino mass, the following topologies are expected: 

• If the gravitino mass is below a few eV/c^, the NLSP decays immediately after its 
production before the tracking subdetectors. In the neutralino NLSP scenario, two very 
energetic photons and missing energy are expected from direct neutralino production. If 
the sleptons are accessible, the cascade decay i ^ ix with two photons and two leptons 
in the final state is a background free search. Nevertheless, two photon searches alone are 
able to exclude most of the neutralino NLSP parameter space in the 'no lifetime' region. 
In the slepton NLSP scenario, if sleptons are directly produced, SUGRA searches for two 
acoplanar^ leptons can be used. Also cascade decays x ^ are possible here, if the 
neutralino pair-production is kinematically allowed. Indeed, its production cross section 
will be greater than the slepton production cross section because of the (3^/ s suppression 
factor on scalars. Finally, in the purely stau NLSP case, where selectrons and smuons 
are heavier than the stau, final states with two leptons and four taus are possible, as 
described in the preceding chapter. 

• For rriQ between a few eV/c^ and a few hundred eV/c^, the NLSP will decay somewhere 
inside the detector (see Fig. |T]^) The decay products of the NLSP, either photons or 
leptons, will present large impact parameters or in the case of slepton NLSP when the 



-'^OPAL and DELPHI also include searches for charginos in e+e~ xtxi ~^ W*xW*x and xtxi ~^ (.viv 
with the NLSP decaying inside the detector But charginos may only be produced at LEP2 if Mm is very 
large. DELPHI has also searched for sgoldstino production and decay to 777 and 17(77, and is able to set limits 
on the scale of SUSY breaking VFq from this search . 

^The acoplanarity a is defined as the angle between the projections of two tracks into the transverse plane. 
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slepton is reconstructed in the tracking volume, the slepton and the lepton will form a 
'kinked' track. 

• For rriQ above a few hundred eV/c^ the NLSP is stable for detector searches as it decays 
outside the detector volume. Only in the slepton NLSP case are direct searches possible: 
highly ionising back-to-back tracks are formed by slow moving sleptons. The situation 
in the neutralino NLSP reverts to the usual SUGRA phenomenology with a neutralino 
LSP escaping the detector, and only searches for indirect production of neutralinos are 
possible. 



4.2 Searches for neutralino NLSP 

Searches for x ~^ 7^ rely on good photon identification. The highly segmented ALEPH 
ECAL makes it possible to measure the energy of isolated photons with an excellent resolution 
(Eq. |2.4|) . The main SM background source for this search is e^e~ uu, from diagrams 
with s-channel Z exchange or t-channel W exchange and one or two initial state radiation 
(ISR) photons. The invariant missing mass M^j^^, = (pe+ + Pe- ~ P'yi "^72)^; will thus have a 
pronounced peak at the Z mass and a long tail from the contribution of the W exchange. 



The simulation of the background was performed using the KK generator |£8[ , which was 
checked with an independent generator NUNGPV Signal processes were simulated using 
SUSYGEN m. 



4.2.1 Acoplanar photons 

Each neutralino is produced with the beam energy A/i/2, and since the gravitino is nearly 
massless and the decay x ~^ 7G is isotropic in the neutralino rest frame, the photon energy 



spectrum is expected to be flat -Emin < E^i, < -Emax, with p2 



E 



I {VS ± - 4m2) (4.1) 



Thus imposing a threshold cut of -Emin = 37GeV on the energy of the least energetic photon is 
very effective against the SM background, where the photons are produced via bremsstrahlung. 
The cut value was optimised for models with a pure bino neutralino and rrie^ = l-lm^ [f70| 



The preselection procedure is based on the requirement of two acoplanar photons and cuts to 
reduce the e"'"e~ — >• 77(7) events, like a < 177° and no additional energy greater than 1 GeV. 
To reduce the doubly radiative Bhabha events e"'"e~ — > e"'"e~77, a cut on the total transverse 
momentum of the photon system pt is very effective. The complete set of cuts for this search 
is summarised in Tab. [4.2|. The distribution of data and background for the missing energy 



and is shown in Fig. 4.1 before the cut on E, 



^72- 



A total of four candidate events are selected in the data while 4.9 are expected from 
background. After full subtraction of the expected background, a 95% C.L. upper limit on 
the neutralino production cross section can be set for an NLSP laboratory lifetime of less than 
3 ns, as shown in Fig. [4.1| . Systematic uncertainties, mainly due to the photon reconstruction 
efficiency and the level of background, only change the cross section limit by less than 1% [^]. 



4.2 Searches for neutralino NLSP 
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10 20 30 40 50 60 70 80 90 50 60 70 80 90 100 

Min(Ej,E2) (GeV) Xi" mass (GeV/c^) 

Figure 4.1: The distributions of (a) the invariant mass of the system recoihng against the photon 
candidates and (b) the energy of the second most energetic photon, for the two photon search before 
the cut on the energy of the second most energetic photon. The SM expectation is the histogram and 
dots with errors are data. On the right, the 95% confidence level upper limit on the e^e~ — > XX cross 
section in the neutralino NLSP scenario when the neutralino lifetime is negligible: < 3 ns. The 
expected cross section for the MGM model [^] (where the neutralino is purely bino and rric^ = l.lm^) 
is also shown. From Ref. [0^ ]. 



4.2.2 Non-pointing photons 

If the neutralino NLSP has medium decay length such that it can decay inside the detector, 
the most probable topology is a single photon which does not point to the interaction region. 
Therefore one of the produced neutralinos is more likely to escape the detector while the other 
decays before the ECAL. Here, the calorimeter granularity is crucial to be able to reconstruct 
the shower axis and thus accurately determine the displaced vertex distance to the IP. The basic 
selection criteria are to select only one identified photon with an impact parameter greater than 
40 cm in the acceptance region | cos^^| < 0.95 and require no charged tracks (which excludes 
photon conversions). To eliminate radiative Bhabha events e^e~ — > e"'"e~7, there must be no 
energy deposited within 14° of the beam axis and less than 1 GeV in additional energy. To 
further reduce ECAL noise and cosmic muons, at least 40% of the photon energy must be 
recorded in the ECAL, there must be no activity in the muon chambers or within a transverse 
distance of 15 cm from the photon candidate in the ECAL. 

The efficiency for this selection reaches a maximum of 10% at neutralino decay lengths of 
around 8 m. After all cuts, 0.8 events are expected from e+e~ — >■ uu'y and 0.2 ±0.2 events from 
cosmic rays and detector noise. Two events are selected in the data. 

A limit on the neutralino mass as a function of its decay length is obtained under the as- 
sumptions of the MGM (minimal gauge mediated) model |7^, with a maximum of 98.8 GeV/c^ 
at 95% C.L. when both two-photon and non-pointing photon searches are included. As will 
be seen in the next chapter, this limit is reduced to ~94 GeV/c^ if the conditions of the MGM 
model are removed. The CDF detector at the Tevatron recorded an event with two opposite 
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Two photons 


Non-pointing photon 


N^=2 or (iV^ > 3 and E^^^^ > OA^/s) 
Accept conversions 
Acoplanarity < 177° 
Ey, > 37GeV 


A^^=l and cIq > 40 cm 
No conversions 
Eu = 

^ECAL > OAE^ and E^^ca? = 


Additional energy < 1 GeV 
cos < 0.95 
Total pt > 0.0375^miss 



Table 4.2: List of cuts for the neutralino NLSP scenario in the no lifetime case (two acoplanar 
photons and ^) and the intermediate lifetime case (non-pointing photon). 



charge electrons with transverse energy of 63 and 36 GeV, two photons of 36 and 30 GeV and 
a total missing transverse energy of 55 GeV [|71|. This type of signature has a very low back- 
ground expectation (~10~'^) and could be interpreted as a GMSB process in the neutralino 
NLSP scenario: qq — > CReR — > e"'"e~xx e"'"e~77GG ||7^. At LEP, the plane {me^,ni^) is di- 



rectly probed by the direct production of neutralinos through the selectron t-channel exchange, 
thus the searches described above are sensitive to the same parameter space. Figure [4.2| shows 
the excluded area in that plane. The region favoured by the CDF event is now completely 
excluded with 95% C.L. 




Figure 4.2: Left, the excluded neutralino mass as a function of the proper lifetime for the two-photon 
search (dark shade) and non-pointing photon search (light shade). On the right, in the neutralino- 
selectron mass plane the 95% C.L. excluded area by photon searches assuming a pure bino neutralino 
NLSP. Overlaid is the region favoured by the CDF event [0], interpreted as qq — > brbr — > e^e~xx ~^ 
e"'"e~77GG [72|. The dark shaded region corresponds to rrie^ < rriy.. From Ref. [34|. 
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4.2.3 Indirect searches 

For neutralinos with decay lengths greater than the detector size, indirect searches have to be 
used to obtain a limit on the neutralino mass. This scenario is identical to gravity mediated 
models with a stable neutralino LSP escaping the detector. Thus searches for sleptons e+e" — *■ 
— > Ix^X ^-iid charginos e+e~ — >• xtXi ~^ W*"'"x W*~X) cis developed for the SUGRA case, 
are utilised here^. The relationship between the x mass and the slepton and chargino masses 
in these processes can be exploited to set limits on the neutralino mass. Searches for two 



acoplanar leptons and missing energy are described in Ref. The exclusion limits in the 



slepton-neutralino mass plane can be seen in Fig. ^2 . 

Charginos in the MSSM can either decay to a W boson and a neutralino (generally referred 
as three-body chargino decay), to a sneutrino and the corresponding lepton (two-body decay), 
or to a slepton and the corresponding neutrino (slepton decay mode). Chargino searches in the 
three-body decay mode are split into three, according to the decays of the virtual W: hadronic 
Xi — > qqx (four jets and missing energy), leptonic Xi ~^ ^^^X (two acoplanar leptons and 
large missing energy) and mixed (two jets and one lepton). The final results are given in 



Ref. WM and the selections are described in Ref 



4.3 Searches for slepton NLSP 
4.3.1 Acoplanar leptons 

If the slepton lifetime is negligible, of the order of a few mm or less, the GMSB signature 
arising from slepton pair-production £ iG does not differ from that of models with gravity 
mediated SUSY breaking and neutralino LSP with a very small mass £ — > • Two energetic 
acoplanar leptons and missing energy are expected. The search used to cover this topology 
is developed in the SUGRA framework but its results can be interpreted in GMSB taking 
~ (the gravitino is at least six orders of magnitude lighter than the neutralino). Details 



of the selection criteria are described in Ref. [^. The largest background is due to W pair 



production followed by leptonic W decay: W ^ iu which is kinematically very similar to 
i — *■ £G for slepton masses around the W mass. Cuts on Mvis and the momentum of the most 
energetic lepton are needed to reduce this type of background^. 

For each lepton flavour, the number of data events selected in year 2000 data is 39, 39 
and 11 for electrons, muons and taus, to be compared with respectively 38.7, 34.7 and 12.2 
expected events from background. The contribution from the WW background is around 80%. 
Thus good agreement between the SM prediction and the experimental results is obtained in 
this topology. The expected and observed mass limits (with WW background subtraction) for 



each slepton flavour are shown in Fig. 4.3 



4.3.2 Large impact parameters and kinks 

Two distinctive topologies with very small background level may arise for intermediate slepton 
lifetimes. If the slepton decays before having produced a reconstructible track, i.e. for A 



^ If the sleptons are lighter than the charginos, two-body decays of charginos will dominate: Xi ~^ ^^'^ or 
even xt ^'^v. 



^ The lepton's energy distribution is flat with the same endpoints i?min,£'max as Eq. 4.1 with m^. Thus 
leptons are expected to be energetic and, specially smuons and staus with only s-channel production, to have 
symmetric distributions in the polar angle 9. 
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Figure 4.3: Excluded regions at 95% CL. in the m^^ versus plane in the acoplanar leptons 
SUGRA search |73]. The observed (shaded area) and expected (solid curve) limits are given for 
BR(£r £x)=l- The dashed curves give the observed limits when slepton cascade decays (into £x2) 
are taken into account, generally for very small neutralino masses. The dashed-dotted curves mark 
the kinematically accessible region in the MSSM at 209 GeV for a neutralino LSP. The table lists the 
lower limits on the slepton masses for ~ 0. For the stau, the minimum limit (shown as the dotted 
curve in the stau plot) is calculated with a mixing angle of 9f ~ 52° which minimises the production 
cross section. 
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Large impact parameter 


Kinks 


Anti-77 


1 or more tracks with 6 TPC hits and pt > 0.03y^ 
^vis > 0.03^5 


Anti-^£, 


Remove conversions 


cosmics, 


ECAL timing with 


single helix fit 


Bhabha 


^vis < 0.65^5 


^vis < 0.90^5 




2 or 4 tracks with > 4 TPC hits 


< 6 tracks with TPC hits. 
Same-charge tracks 


Tracks 




with dr(i, < 5 cm and d^ < 20 cm. 


quality 


Track ^1: do > 1cm and pi > 0.01-y/s 


Inner track: d^j^ < 0.5 cm. 


cuts 


Track #2: d'^ > 0.025 cm 


(2) 

Outer track: dg > 1 cm 
with at least 4 TPC hits, 
p > 0.015^5 and a2,z > 18.2° 


Anti-77, 


Acoplanarity < 175° 


ai,2 > 5° or ai,2 > 10° 


U and 




if outer track has no TPC hits. 


had. ints. 


AcoUinearity < 11.5° 


El < 5 GeV 



Table 4.3: List of cuts for the large impact parameter and kinks selections. 



between a few mm and a few cm, the lepton tracks will present large distances to the beam axis 
(impact parameters). Otherwise, if the slepton reaches the tracking devices and is reconstructed 
as a charged track decaying with A > 40 cm, the typical signature consists then of kinks formed 
between the slepton track and its lepton decay product. 

The preselection procedure is the same for both large impact parameter and kinks searches. 



Anti-77, dilepton and cosmic rays cuts are applied as listed in Tab. 4.3 



Search for large impact parameters 

Events containing exactly two or four charged tracks are selected to attempt tau identification. 
At least one track (track 1) is required to have an impact parameter greater than 1 cm and 
total momentum of at least l%^/s. Only one or three additional tracks with at least four 
TPC hits are allowed. If there are three tracks, they must be consistent with a three-prong 
tau decay. These tau triplets are treated as one track (track 2) calculating the mean value of 
the three impact parameters as the tau dQ and the sum of the momenta as the tau direction. 
Track 2 is required to have an impact parameter of at least 0.025 cm. The acoplanarity and 
acollinearityQ of tracks 1 and 2 are used to suppress dilepton and 77 background. 



Search for kinks 

Kinks are reconstructed by searching for same-sign tracks crossing in the r0 plane (transverse 
to the beam axis) or approaching each other in this projection closer than 5 cm inside the 
tracking volume. At this point, the separation between the two in the z direction should be 
less than 20 cm. 

In the search for kinks the most important sources of background are hard bremsstrahlung 
and hadronic interactions in the material of the detector from Kg, K^ or vr^ decays. To reduce 



^The acollinearity is defined as the opening angle between two tracks. 
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the former, the angle between the primary and secondary track ai^2 must be greater than 
5° at the kink vertex: kinks from bremsstrahlung are peaked at low values of 0:1.2. Finally, 
the residual background comes from ditau events with hadronic interactions in the material, 
i.e. when the kink is found between 26 and 34 cm radius. These events are characterised by 
accompanying coUinear particles to the interacting hadron, thus the energy deposited in the 
calorimeters within 5° around the axis from the IP to the kink must be less than 5 GeV to 
reject these events. 

Combination 

Combined efficiencies above 20% and with a maximum of ~65% for stau samples are obtained 
in the NLSP lifetime range covered by the large impact parameter and kinks searches from 
~10^^'^ to 10^^ s. The efficiency is higher for selectron and smuon samples due to the absence of 
three-prong decays and the higher momenta of the leptons. In the 189 — 209 GeV data sample 
only one event is selected by the combined search with 1.1 expected from SM backgrounds. 
The event is recorded at 189 GeV and is compatible with a three-prong tau decay in the large 
impact parameter selection [ [f6[] . 

4.3.3 Heavy stable charged particles 

Long-lived sleptons can be detected as two back-to-back charged particles with unusually high 
values of dE/dx for their momentum. They are not expected to interact hadronically (they 
do not have colour charge) nor electromagnetically (since they are very heavy). Thus they 
should be very similar to dimuon events only distinctive by their high specific ionisation. 
In the kinematical limit, where the sleptons are expected to be very slow, the energy loss 
becomes high enough to saturate the TPC electronics. To cover this interesting mass region 
an additional selection based on calorimeter information and the pattern of saturated hits was 
developed . 



Stable sleptons 



iVch = 2 each with > O.l^bcam 
\cos6\ < 0.9 with more than 1 ITC hit 
Moil + M02I < 0.3cm 
\zoi \ + \zo2\ < 5cm 
No identified electron 
Acollinearity > 160° 

bi -P2I < 3W(T, 



pi 



0", 



p2 



< 250 MeV 



-Eecali + -Eecal2 < 20 GeV and 

^HCALl + ^HCAL2 < 50 GeV 



0.52 < mi,2/^beam < 0.80 

Acollinearity > 174° 



0.80 < mi,2/^beam < 0.98 

Rii + Ri2 > 10 



Table 4.4: List of cuts for the heavy stable charged particle selection for intermediate-high masses. 
The region close to the kinematic hmit (mi^2 > 0.95£'bcam) is explored by the search for saturated 



TPC hits described in Ref. [^|. 
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Figure 4.4: (a) Efficiency of the search for pair-production of massive, spin-0 stable particles as a 
function of m/-Ebeam with (circles) and without (squares) analysis of events with no good tracks due to 
saturation, (b) The maximum allowed cross section calculated at 95% confidence level with data taken 
in year 2000 at centre-of-mass energies of 202-209 GeV (filled area) with the heavy stable charged 
particles search (including the analysis of saturated hits). The solid lines represent the minimum 
expected cross section times branching ratio for the lightest and heaviest f pair-production at ^/s = 
207 GeV. Thus stau masses of up to 97GeV/c^ are excluded by this search alone in the stau NLSP 
scenario. 



Table describes the cuts used in the selection of medium-high masses. The track quality 
criteria, the first four cuts in Tab. |4.4| , are imposed to reject events with poorly reconstructed 
tracks and cosmic rays. The electron ID veto and the angle acceptance reduce the large Bhabha 
background. The acollinearity cut effectively reduces the radiative Z returns, 77 and t^t~ 
events. The momentum cuts also reject 77 and rr events. Requiring no photons with energy 
greater than 250 MeV considerably reduces on-shell Z events and Bhabhas. Finally the ECAL 
and HCAL energy deposits help in Bhabha rejection. After these cuts are applied, the high 
mass of the signal has not been used and the background is dominated by events. If the 

reconstructed mass of the particles = \/ -^beam ~ Pi momenta of the tracks are 

high enough to produce an energy loss similar to that of ordinary particles. Only in the high 
mass range will the dE/dx information be useful in terms of the estimator Rj = {I — (/^))/o"/. 
Here one compares the measured dE/dx, I, to that expected for a muon (J^), where aj is the 
expected resolution of the measurement. 

The case of masses rrij above 0.95i?beam is described in Ref. where the selection makes 
use of the position of the two most energetic calorimeter objects and the IP to perform a 
single-helix fit of the transverse momentum and polar angle of the two track candidates. Thus 
the trajectories of the particles can be estimated without ITC or TPC information and are 
then compared with the saturated and unsaturated hits. The efficiency ranges roughly from 
50 to 70% for the intermediate-high mass, between 0.5 and 0.93 m/_Ebeam- If the analysis of 
saturated hits is included as described in Ref. |7^, the efficiency is improved and maintained 
with an approximately flat distribution (within ±10%) reaching the highest values of m/i?beam- 
The efficiencies for both medium-high and ultra-high masses selections are shown in Fig. [4.4| a. 
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Figure 4.5: Excluded stau mass at 95% confidence level as a function of its lifetime (shaded area) 
from direct searches. Dashed curves give the limits from the different topologies. The search for 
acoplanar leptons covers the case of small lifetimes, searches for tracks with large impact parameter 
and for kinks are used in the intermediate range, whereas for large lifetimes a search for heavy stable 
charged particles is performed. The dotted curve gives the expected limit. From Ref. ||65[| . 



The combination of intermediate-high mass searches and the saturated hits analysis selects 
one data event while 2.3 are expected from background, mainly dimuons. Systematic effects 
have been studied with independently selected dimuon samples from the data to check the 
performance of the tracking system. Specific ionisation effects have been checked with electrons, 
muons and pions over a large range in momentum. A total systematic error of less than 5% 



is estimated and conservatively applied to the selection efficiency. Figure |4.4| b shows that 
cross sections as low as 0.02 pb are excluded using only the year 2000 data set for stau pair- 
production. 

Combination of direct slepton NLSP production searches 

When all the above selections are combined including all lifetimes, the 95% confidence level 
lower limits on the right-slepton masses, independent of lifetime, are set at 83, 88 and 77GeV/c^ 
for selectron, smuon and stau, respectively. The selectron mass limit is obtained neglecting the 



t-channel exchange contribution to the production cross section. Figure [4.5| gives the excluded 
stau mass as a function of its lifetime when all the above described searches are combined. The 
dotted lines represent the independent limits from searches for negligible lifetime (acoplanar 
leptons from SUGRA with ~ 0), intermediate lifetime (large impact parameters and kinks) 
and long lifetime (heavy stable charged particles). The .OR. of the searches is the solid line [|65|] . 
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4.3.4 Cascade decays of neutralinos 

Searches for slepton NLSP pairs when they are produced directly at LEP and decay to a 
lepton and a gravitino have been described above. However, if the XX production is allowed, 
it generally provides the dominant signal in the slepton NLSP scenario. This is because the 
production cross section for neutralinos is larger than for sleptons since it is not suppressed, 
even for > (see Fig. |3.2| b). The expected signal is therefore four leptons and missing 
energy from the cascade decay: e"'"e~ — > XX ~^ £'£'G. Two of the leptons may be 

very soft depending on the mass difference between the neutralino and the slepton. However, 
in half of the cases the two hard leptons may have the same charge because of the Majorana 
nature of the neutralino. This makes this type of decay a very clean discovery signal in models 
with small N5 and tan j3 not too large. In ALEPH, searches for ^ have been carried out 

for both negligible slepton lifetimes and intermediate slepton lifetimes, when two of the tracks 
present kinks or large impact parameters Since the neutralinos decay independently of 



one another, there are six possible topologies for the slepton co-NLSP, namely: ee, /i/x, rr, e/x, 
ef and jlf. In the stau NLSP case, only the ff topology is relevant. 



Selection 


ee /i/i TT e/i er /ir 


Negligible Lifetime 


obs. 
exp. 


6 1 22 2 8 5 
5.15 0.44 16.49 3.65 5.52 5.72 


Intermediate Lifetime 


Short A 


obs. 
exp. 


5 

5.25 


Long A 


obs. 
exp. 


4 
1.51 



Table 4.5: Number of observed and expected events in the searches for four leptons and missing 
energy in the 189 — 209 GeV data. In the intermediate lifetime case there is large overlap between 
the different selections and the results are given for the total. 



Four leptons and ^ 

In the 'no lifetime' case for topologies not involving taus, the presence of at least three identified 
electrons or muons is required but not more than four charged tracks. This helps to reject 
qq and four-fermion events with high multiplicity. For each topology, at least two energetic 
leptons are required which must be acoplanar, acollinear and their energies must be in the 
range allowed by the signal kinematical properties. These cuts reject effectively the dilepton 
background. The kinematic cuts reduce also leptonic WW decays and Wez/ background. 

For topologies involving taus ef and /if, the number of charged tracks per event and the 
mean number of neutrinos carrying away missing energy and momentum are considerably 
bigger than in ee, jljl and e/i topologies. Thus four, five or six charged tracks are allowed with 
at least one identified lepton. Finally, for f f topologies at least three identified tau jets are 
required with not more than two other charged tracks. Cuts on the Durham cluster thresholds, 
the energy of the event and the momentum of the most energetic lepton are used to reduce 
the main background contributions, WW and rr events. The selection cuts are described in 
detail in Ref. [|65| . 

Efficiencies of up to 80% are achieved in the ee and fijl selections. For topologies with 
taus, the efficiency can be as low as 30% for small differences between the neutralino and the 
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stau, reaching a maximum of 60% for energetic taus. Figure [4.6^ displays the efficiencies as 
a function of Am^^. The number of selected events in the 189 — 209 GeV data is shown in 



Tab. 4.5 for each topology. 



Four leptons and ^ with lifetime 

A detailed description of this search can be found in Ref. |7^. The topology consists of two 
tracks with large impact parameter or kinks and two tracks originating from the interaction 
point. Thus the basic selection requirements rely on tracks with more than four hits in the 
TPC and high-o?o- The main background sources are nuclear interactions, photon conversions, 



ECAL 'splashbacks' and cosmic muons. Table [4.5| shows the comparison between expected 
background and selected data events in the large impact parameter and kinks searches. All 
bar one of the selected events can be assigned to a known source of background. 

The efficiency of the long decay selection (kinks) as a function of lifetime is presented in 
Fig. for the ff signal. A peak of 80% efficiency is achieved at decay lengths around 
10 cm deteriorating for small Am^f = — rrif. Also shown in that figure is the efficiency for 
the negligible lifetime selection when applied to signal data with slepton lifetime. Overall, an 
efficiency in excess of 10% is maintained for slepton decay lengths from ~1 mm to ~10m for 



all channels [78 . 




5 10 15 20 25 10"^ 1 10 10^ 10^ 

Am (GeV) Slepton decay length (cm) 



Figure 4.6: (a) Selection efficiencies for the six different event topologies versus ^Tn^i = ~ 
in the negligible slepton lifetime case. The spread of points observed for a given topology is due to 
different values of neutralino and slepton masses, (b) Probability for a signal stau-pair event to be 
selected by at least one of the six topological searches versus slepton decay length. The set of curves 
with higher efficiency in the 0.1cm area corresponds to the prompt decay selection. Those peaking 
at ~10cm correspond to the short and long decay length selection. Different lines correspond to 
different points in the (m^,mf ) space as indicated. 



4.4 Neutral Higgs boson searches in the MSSM 
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4.4 Neutral Higgs boson searches in the MSSM 

Searches for neutral MSSM Higgs bosons can also be used to constrain the GMSB parameter 
space. Although these searches are performed in the framework of gravity mediated scenarios, 
the Higgs sector is intrinsically the same in GMSB models. Thus limits from one can be applied 
to the other taking into account the different parameters involved in the determination of the 
Higgs masses and couplings. 

At LEP2 the lightest CP-even and CP-odd neutral Higgs bosons can be produced 
through Higgsstmhlung ^ Ta* ^ KL and associated production e^e~ ^ Z* ^ hA (see 

Fig. |1.6| d). The cross sections of these two reactions are complementary, Higgsstrahlung is 
dominant at low tan/3; while hA has larger cross sections at large tan/?. This is due to the 
MSSM suppression factors sin^(/3 — a) and cos^(/3 — a) |]^: 

cr(e+e- hZ) ~ sin^(/5 - a)asM (4.2) 
a{e^e~ hA) ~ cos^(/3 — Q;)crsM (4.3) 

where a^u is the production cross section for the SM Higgs boson radiated off from a Z, and 
a is the mixing angle of the two Higgs fields that produce h and H . 

From now on h denotes the MSSM as well as the SM Higgs boson. The two most relevant 
decays of h and A are to hb and r+r~. Thus four main signal topologies arise: four jets, missing 
energy, leptonic and tau final states, as listed in Tab. |0. 



Topology 


Z decay h decay A decay 


BR(/iZ^ ) BR{hA ) 


Four jets 
Missing energy 
Leptonic 

With taus 


qq hh bb 
vv bb XX, GG, . . . 

/+/- bb 
qq TT bb 

TT bb TT 


70% 86% 
20% small 
7% 

3% 13% 



Table 4.6: The four main topologies considered in neutral Higgs boson searches, as derived from 
hZ or hA decays. The branching ratios for hZ decays are given for the Standard Model Higgs with 
mass 115GeV/c^. In hA decays, the hsted branching fractions are valid for maximal h masses and 
tan f3 < 30. 



Searches for all the above topologies have been performed in ALEPH and the selections are 
described in Ref. []79[ with the final results given in Ref. here only a general description 
is given. 

The four-jet final state is the most likely signal given its large branching fraction, but suffers 
from irreducible backgrounds such as ZZ, bb with FSR hard gluon splitting to a pair of 6-quarks 
and WW with 6-quark misidentification. Nevertheless, the four jets in the signal are expected 
to be contained in a plane if the Z and h or h and A are produced almost at rest. In the first 
case, the invariant mass of two of the jets should be close to mz and the other two contain b 
flavour, while for the second case all four jets contain B hadrons. 

The missing energy final state is characterised in hZ production by Mmiss ~ and sub- 
stantial missing transverse momentum. In the case of hA production, if other SUSY particles 
are kinematically accessible or the bosons undergo a WW fusion process into bbv^v^, a similar 
event topology is obtained. 
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If the Z decays to electrons or muons, the two leptons must reconstruct the Z mass and the 
two jets should be tagged as 6-jets. mh can then be reconstructed accurately from the system 
recoiling against the two leptons with a typical mass resolution of about 1.5 GeV/c^. 

The relevant limit for GMSB searches is obtained by the combination of the /iZ and the 
hA searches in the plane given by the h mass and the factor sin^(/5 — a). This plane is model 
independent and the present ALEPH exclusion is shown in Fig. [4.7| . 




Figure 4.7: The 95% confidence level exclusion contours for the /iZ and hA searches as a function of 
siv?{j3 — a) (dashed lines). The combined exclusion is shown by the hatched area and the dotted line 
indicates the expected exclusion. From Ref. l80|] . 



Chapter 5 

Interpretation of results 



The previous two chapters describe all the searches performed in the ALEPH collab- 
oration in the search for GMSB topologies. No hint for a signal is observed in any of 
them. This chapter is therefore devoted to understand what the lack of evidence for 
such signals may tell about the theoretical models that predict them. Understanding 
the effect of 'not seeing' any signal at LEP is important to design new searches at 
future colliders and maybe address new theoretical models. 

The interpretation of the search results in this work has been limited so far to a 
specific model or constrained by simplifying assumptions made to derive the limits. 
To extrapolate these results to the more general case, a scan is performed over the 
six parameters that fully determine GMSB phenomenology. It is then possible to 
test the robustness of the previous results, checking how they may change when 
no assumptions whatsoever are made, and to interpret them in a wider context. 
Moreover, each of the above searches will contribute to exclude different regions in 
the minimal GMSB parameter space, thus constraining the different parameters of 
the model distinctively. It will be important to determine which searches set the most 
effective limits on the theoretical parameters. 

The final ALEPH combined limits on the parameters set in this chapter are in most 
cases very close to the sensitivity of the detector to such signals. A discussion on the 
allowed ranges of the parameters is presented, to asses the validity of these limits if 
the scan ranges were enlarged. 
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5.1 The scan 

The entire sparticle mass spectrum and all possible phenomenological topologies can be deter- 



mined by only six parameters in GMSB models |32|. A scan over these parameters has been 
performed to interpret the results of the experimental searches. The aim is to understand 
which searches contribute to exclude which regions in the parameter space and to be able to 
set a lower limit on the mass of the NLSP and on the universal mass scale A, independent 
of the NLSP lifetime (i.e., for all gravitino masses). The program ISASUSY 7.51 |^ was 



used to calculate masses, cross sections, and branching ratios for all SUSY particles at each 
point in the parameter space. Radiative corrections to chargino and neutralino masses were 
applied. The Higgs bosons' masses and couplings were calculated using the improved Haber 
et al. treatment at two loops implemented in HZHA 3.0 ||5T| with a top mass of 175GeV/c^. 



5.1.1 Scan ranges 



The ranges of the scan are given in Tab. for each parameter as well as the step size. In 
total, over 2.3 million points in the minimal GMSB parameter space have been tested. If the 
NLSP mass is greater than the highest beam energy (104.5 GeV) the point is not considered 
(and cannot be excluded). The scan is not exhaustive but it covers a large portion of the 
allowed range of the minimal GMSB parameter space, as discussed next. 



Parameter 


Lower limit Upper limit Steps 


M,n 

rriQ 
A 

N5 
tan/5 
sign(/i) 


10^ GeV/c2 10^2 GeV/c2 7 
10-1 eV/c^ 10^eV/c2 15 
10^ GeV/c^ min(v^, M^) 100 
1 5 5 
1.5 40 22 
+ 2 



Table 5.1: Minimal set of parameters and their ranges of variation used in the scan. 



Validity of the ranges 

Mm An approximate upper bound on the common messenger mass is imposed by flavour 
universality. This can only be maintained if gravity mediated contributions to the ob- 
servable masses (of the order Fo/Mp) are kept small ^O{10~^). Thus to avoid large 
flavour violating transitions 



1 (T 

Mn, < — — Mp ~ 10^^ GeV/c' (5.1) 
IO2 47r 

Furthermore, NLSP decays in the early Universe must not spoil the predictions from 
Big Bang nucleosynthesis. These bounds can further reduce the upper limit on M^ as 
described in Ref. |]5^. The lower limit is determined by A with the condition M^ > A. 



ttIq The gravitino mass range is deduced from Cosmology and from indirect searches. Direct 
collider searches e"'"e- — > Gx set a lower limit on ttiq of ~10-^eV/c^ (see Sec. |1.5| ), but 



5.1 The scan 
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previous indirect limits suggest that > 10 ^eV/c^ In any case, the gravitino 



mass controls the lifetime of the NLSP. Phenomenologically, setting ttt-q below 10~^ eV/c^ 
would not improve the scan, since that is enough to cover all cases with negligible NLSP 
lifetime ~10~^^ s. As regards the upper limit of lOOkeV/c^, it is imposed by standard Big 



Bang constraints [£3|. A heavier gravitino would overdose the Universe at an early stage, 
because p3/2(T) = mQUQiT) > Pc(T) where is the number density of gravitinos, which 
depends on scattering processes off thermal radiation and the decay of heavier particles 
into gravitinos, and pc is the critical density. Even assuming dilution mechanisms such 
as inflation, an overproduction of heavier gravitinos in the reheating epoch can lead to 
cosmologically unacceptable values of the relic density 

A The universal mass scale parameter is constrained from above by the mass of scalar mes- 
sengers. 

{mm^ = ^ > ^ < (5 3) 

(5.3) 



This upper limit A < min(Mni, a/FV) means that values up to lOOTeV/c^ were scanned. 
Higher values would drive the NLSP mass beyond 1 TeV/c^. 

N5 The number of messenger families is allowed to vary up to five. The messenger sector 
affects how the gauge couplings unify at the GUT scale. In order to maintain corrections 
to the GUT gauge coupling small, and thus allow for unification, N5 and should 




obey |2|: 

As can be seen in Fig. pT7[b, for Mm as low as lOOTeV/c^ N5 can be at most five. Yet, 
for Mjn = 10^°GeV/c^ N5 as large as ten is allowed. However, values of N5 = 6 do not 
change significantly the results as will be proved later. 

tan (3 The ratio of the two Higgs doublets vacuum expectation values has been changed from 
1.5 up to 40 in the scan. Lower values of tan/? are strongly disfavoured by Higgs boson 
searches, while calculations of SUSY parameters in the ultra-high regime, beyond ~40, 
become unreliable. 

sign(/Lt) The value of the mass mixing parameter between the two Higgs doublets can be 
computed by enforcing correct electroweak symmetry breaking (Eq. |1.14|) , thus only an 
ambiguity in its sign remains. 



5.1.2 Exclusion procedure 

Once the possible GMSB models are available, the exclusion procedure takes place: 

• First, the contribution from unknown processes to the Z width cannot exceed 4.7MeV. 
The extremely precise measurement of Fz from LEPl allows one to exclude every new 
particle that could contribute to the Z decay width beyond the small experimental error. 
The procedure to determine AFz is sketched in Fig. ^.1| . The Z width depends on the 
Higgs mass m/j, the strong coupling constant a^, the electroweak coupling constant 
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SM lower bound-. 



-^95:5 



Figure 5.1: Procedure to determine the allowed SUSY contribution to Tz- 



and the top- and bottom-quark masses. Thus all these parameters, and their experi- 
mental errors, enter in the calculation of the predicted Fl"^ ~ 2496.1 ± 1.2 MeV, which 
defines a one-sigma lower bound (2494.9 MeV). This was calculated with the program 
ZFITTER and the latest experimental values of the parameters [^]. Then a Gaussian 
distribution is assumed centered on the experimental value F^^^ ~ 2495.2 ±2.3 MeV. The 
95% C.L. upper bound is calculated by dividing the remaining area under the gaussian 
curve in the ratio 95:5, to obtain r^^^^^ = 2499.6 MeV. The difi'erence between the SM 
lower bound and the 95% C.L. upper limit r^^''^^ is AFz = 4.7 MeV. 

• Then the LEPl limits are applied. The upper limit on the production cross section 
derived from heavy stable charged particles ergs = 0.296 pb |SB| was used. Which generally 
means that NLSP masses below ~mz/2 are excluded by LEPl. 

• Finally, cross section limits from GMSB searches (and SUGRA searches for sleptons and 
charginos) as described in the preceding chapters are imposed on every point for which 
the relevant topology is accessible. In addition, limits from neutral Higgs boson searches 
are imposed on every scan point taking the exclusion in the plane (sin^(/? — a),mh) from 
Fig. 

Thus a point in the parameter space is excluded if it is kinematically accessible and the 
excluded cross section from at least one search is less than the production cross section for 
that topology, taking into account the branching ratio. Cross sections for SUSY particles were 
calculated at seven different averaged centre-of-mass energies, as listed in Tab. PTT| . 

5.2 Lower limit on the NLSP mass 



In the stau NLSP scenario, the 95% C.L. limit on ruf set by direct searches (Fig. [4.5|) is almost 
unaffected by the stau mixing and is 77GeV/c^ for any stau lifetime, over the full scan range. 
Thus Fig. |4.5| from direct slepton production searches sets the minimum allowed mass for the 
lightest stau as a function of the stau lifetime, over the full scan range. Searches in the stau 
NLSP scenario for neutralino production (for four leptons and missing energy, with prompt or 
short decays) are not able to improve this limit. 



5.2 Lower limit on the NLSP mass 



87 




log.JT /S) 



Figure 5.2: Absolute lower limit on the neutralino mass in the x NLSP scenario from GMSB searches 
as derived from the scan. The dashed line gives the lower limit obtained when searches for /iZ and 
hA decays are included. The areas covered by searches for photons, and sleptons and charginos in 
SUGRA, are also shown. The search for single non-pointing photons extends into very long NLSP 
lifetimes, since the probability for only one neutralino decaying inside the detector is still large for 
lifetimes ~300m. 



Limits on in the neutralino NLSP scenario were shown in Fig. [4.2| assuming mg = l.lm^. 
When this condition is relaxed and the full scan is tested, < 94GeV/c^ remains excluded 
for short neutralino lifetimes. The dependence of this limit on the neutralino lifetime can be 
seen in Fig. 5^ where searches for sleptons and charginos have been included to cover the 
long lifetime case. The absolute limit independent of lifetime is set at 54 GeV/c^ by indirect 
searches for sleptons and charginos in the long lifetime case. 

The interplay of the different searches in the (m-^, irtf) plane is shown in Fig. |5.3| . For short 
NLSP lifetimes (Fig. |5.3| a) searches for two photons, two leptons, four leptons and six leptons 
contribute to exclude points in this plane. It can be seen here how searches for four leptons 
and missing energy are able to extend the limit from acoplanar leptons searches in the slepton 
NLSP scenario up to 94GeV/c^ in neutralino mass. Nevertheless, the search for six leptons 
and missing energy, described in Chapter |^ in this work, covers most of that area and is also 
able to exclude fi masses below 84 GeV for < 130 GeV. 

In Fig. 5^b the case of long NLSP lifetimes is presented. Limits in the neutralino NLSP 
region are less constraining than those for the short neutralino lifetimes, due to the use of 



88 



Interpretation of results 



ALEPH Vs= 1 89-209 GeV ALEPH Vs= 1 89-209 GeV 




2 2 

X mass (GeV/c ) % mass (GeV/c ) 



Figure 5.3: Region excluded by tlie different analyses described in the text at 95% confidence level 
in the (m^,mf) plane for (a) short NLSP lifetimes (m^ < lOeV/c^) and (b) long NLSP lifetimes 
(mg > IkeV/c^). Points in the dark region are not accessible to the scan. The absolute NLSP mass 
limit is set at 54GeV/c^ in b by the intersection of chargino and slepton searches. 



topologies with indirect neutralino production. The absolute lower limit on the NLSP mass of 

^NLSP > 54GeV/c2 I (5.5) 



is visible in this plot, determined by the chargino and sleptons searches. This point is found 
at Ns = 1, tan^g = 3, A = 39TeV/c2, M^, = 10^° GeV/c^ and = lO^eV/c^, where the 
neutralino is the NLSP with the i masses around 96GeV/c^ and all other supersymmetric 
particles above threshold. 

The impact of the neutral Higgs boson searches on the neutralino and stau mass limits 
is shown in Fig. as a function of tan /3. In this case, the NLSP absolute mass limit is 
77GeV/c^ obtained for large tan/? and in the stau NLSP scenario. Small values of tan/3 are 
excluded almost independently of N5. 



5.3 Lower limit on the mass scale parameter A 

The parameter A represents the energy scale at which the messenger particles couple to the 
visible sector and hence fixes the universal mass scale of SUSY particles. Hence a lower limit 
on the NLSP mass implies a lower limit on A. Since the relation between the masses and A 
depends on N5 (see Eqs. [1.17| and |1.18|) , it is interesting to see how this limit changes with 



N5. The excluded values for the parameter A as a function of tan/3 are shown in Fig. |5.5| for 
different values of N5. In the short lifetime case the lower limit on A appears around 12 TeV/ 



5.3 Lower limit on the mass scale parameter A 
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Figure 5.4: Lower limits on the masses of (a) x and (b) fi as a function of tan/3, for different values 
of N5, as set by the Higgs boson searches. The shaded area represents the minimum excluded area, 
for any N5, as derived from GMSB searches alone. 



for N5 = 5. The absolute limit from the full scan is set at 



A > lOTeV/c^ 



(5.6) 



in the long NLSP lifetime case. This limit is set at N5 = 5, tan (3 = 1.5, and Mm = 10^^ GeV/c^. 
The neutralino is the NLSP here with a mass of 73GeV/c^, slepton masses are around 
76GeV/c^ and all other sparticles are above threshold. The excluded A increases with tan/3 
as slepton NLSP searches become more relevant and improve on the indirect limits from neu- 
tralino searches alone. 



Figure [5.5| c shows the lower limit on A for any NLSP lifetime as a function of tan When 
searches for neutral Higgs bosons are included, the excluded area significantly increases, as is 
shown in Fig. |5T5|d, specially for low tan/?. For example, for N5 = 1, A up to 67TeV/c^ and 
tan /3 up to 6 are excluded. 

Finally, the lower limit on A independent of lifetime and tan/5 is shown in Fig. ^a as a 
function of N5. If the Higgs exclusion is taken into account the limit on A increases up to 
66, 39, 26, 20 and 16 for N5 = 1, 2, 3, 4, 5 respectively, assuming = 175 GeV/c^. 

The implication of possible extensions of the scan ranges is discussed next. 



5.3.1 Validity of the limits 

Two additional small scans are performed. One to check the limits if the number of families is 
increased to six, and the other to see the effect of a 180 GeV/c^ top mass on the limits derived 
with mt = 175GeV/c2. 

The lower limit on A for N5 = 6 is 9.3 and 14.4TeV/c^ for GMSB searches alone and 
Higgs boson searches, respectively. Thus the limit on the mass parameter A is reduced if more 
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tanp tanP 




tanP tanP 

Figure 5.5: Region excluded at 95% C.L. in tlie (A,tan/3) plane for (a) siiort, (b) long and (c) any 
NLSP lifetime. The impact of the Higgs search is included in d. Values of tan (3 less than 3 are 
excluded for large N5 at any mass parameter A, while tan/? up to 6 can be excluded for N5 = 1. 



5.4 Lower limit on the gravitino mass 
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families are included in the messenger sector. But it does so by a small amount, from around 
10 to around 9TeV/c^, compared to the decrease from four to five messenger families. 

The top mass is measured directly to be 174.3 ± 5.1GeV/c^ [Q. The general scan was 
performed assuming 175 GeV/c^. If this is increased by the experimental error to 180 GeV/c^, 
the lower limit on A for N5 = 5 from neutral Higgs boson searches is reduced from 16 to 
15TeV/c2. 

A change in the ranges of other parameters would not affect these limits on A because they 
are independent of the NLSP lifetime (mg is completely covered) and larger values of tan/? 
(> 40) may only produce lighter fi which are already excluded. 



5.4 Lower limit on the gravitino mass 



The equation that relates the gravitino mass to the scale of SUSY breaking Fq can be exploited 
to put an indirect limit on the gravitino mass. Using Eq. and the fact that the messengers 
'feel' the SUSY breaking at a lower scale: Fq > Fm, 



Fo 



A2 



(5.7) 



Thus the lower limit on A can then be converted into an indirect limit on the gravitino mass. 
The dependence of rriQ on N5 is illustrated in Fig. p.6| b, and the lower limit of A > lOTeV/c^ 
(16 TeV/c^ if the Higgs limits are included) implies a lower limit on rriQ of 0.024 (0.061) eV/c^. 
Therefore, 

' > 2 X 10-2 eV/c2 I (5.8) 



This limit is set for N5 = 5, but for only one messenger family and, if the hA and hZ limits 
are included, it is increased to leV/c^. 
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Figure 5.6: Exclusions at 95% confidence level (a) for A and (b) for ttIq as a function of N^, derived 
from the minimal GMSB scan (shaded). The unshaded bars represent the excluded region when the 
neutral Higgs boson exclusion is applied. 



Chapter 6 

Summary and discussion 



6.1 This work 

Models of super symmetry where the MSSM masses are generated by radiative corrections 
due to messenger particles coupling through the gauge interactions offer an extremely rich 
phenomenology. Several differences arise with respect to the usual approach of gravity mediated 



SUSY breaking, as listed in Tab. 



mSUGRA 


mGMSB 


mi/2, mo, Ao, tan/3, sign(/i) 


Mm, m^, A, Ns, tan/3, sign(/i) 


Gaugino unification mass relations (Eq. |1.16|) 
Radiative EWSB (Eq. |1.14|) 


~ 10^^ GeV 
Mm ~ Mp 
possible FCNC 
Aq = free parameter 
ml ^ ^1/2 ^ — > 

ml <^ m^y2 ^ > TTT'l — l^x 

~ 100 GeV/c2 
X? LSP 


~ 10^ - 10^° GeV 
Mm ^ 10^ - 10^5 GeV 
no FCNC 
v4 ~ 

niq : m^^ : m^^ : rriy. = 11 : 2.5 : 1.1 : V'Ns 

TTT-Q < IkeV 
G LSP with xi or ii NLSP, 
NLSP has lifetime 



Table 6.1: Comparison between minimal supergravity and minimal gauge mediated models. 



If SUSY breaking is generated at Mm as opposed to Mp, the gravitino mass is drastically 
reduced and it becomes the LSP. The decay width of the next lightest sparticle to the grav- 
itino Pnlsp oc m^i^gpm'^^Mp'^ can then be small enough to produce displaced vertices inside 
detectors. The nature of the NLSP can be read off from Tab. |6.1| to be the neutralino if only 
one family of messengers exists, or the lightest slepton (usually the stau) in any other case. 

In this last type of scenario, where the fi is the NLSP and the other sleptons may be even 
heavier than the lightest neutralino x, signatures with two soft electrons or muons, four taus 
and missing energy are expected in ALEPH. A search for this topology has been developed 
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and, for the first time, special attention has been placed on three-body decays of selectrons and 
smuons to the corresponding lepton, a tau and the stau NLSP. A dedicated selection covered 
the region of very soft leptons, when the mass splitting between the sleptons and the NLSP 
is small. This search alone conservatively excludes selectron and smuon masses of up to 94 
and 95GeV/c^ respectively, with 95% confidence level in the case that the stau lifetime is 
negligible. The existing limits for selectron and smuon NLSP direct production were 98 and 
93GeV/c^, respectively. 

But the importance of this search becomes apparent when a full scan of all possible GMSB 
models is performed and the very many different searches are combined to understand the 
exclusions in a wide context. It is then proven to significantly exclude an as yet underestimated 
area in parameter space over that ruled out by cascade decays of neutralinos to sleptons and 
direct slepton decays (see Fig. ^73|a). Searches to cover almost every possible signature have 
been developed in ALEPH and are detailed here, including the neutral Higgs bosons searches. 
The interplay between the parameters of the models and the sensitivity of the limits to its 
variation has been studied. From the combined exclusion, and allowing the models to extend 
in the parameter space as much as physically possible, three model independent limits are set: 

• The mass of the NLSP is excluded up to 77GeV/c^ (see Tab. for details). 



NLSP 


mass limit (95% C.L.) 


validity 


X 


92 GeV/c2 
54 GeV/c^ 


short X lifetime {rriQ < lOeV/c^) 
any lifetime 




77 GeV/c2 


any lifetime 


any 


77 GeV/c2 


Higgs exclusion 



Table 6.2: NLSP mass limits from GMSB searches, as derived from the scan. 



• The universal mass scale parameter A is excluded up to 16 TeV/c^ for N5 = 5. This limit 
is not expected to be affected much for larger N5. Values of A in excess of ~ ^ TeV/c^ 
would produce very heavy sparticles, destroying the solution to the hierarchy problem. 

• The mass of the gravitino can be indirectly excluded up to 6 x 10^^ eV/c^, or equivalently, 
the SUSY breaking scale \/F^ has to be greater than 16TeV. Cosmology bounds v^F^ 
from above at ~1000TeV. 

The above limits are calculated with neutral Higgs bosons results and GMSB results as 
described in Chapter |[ 



6.2 GMSB after LEP 



Supersymmetry has escaped detection at LEP, but the techniques and limits from this machine 
will be valuable for future experiments which will surely discover it, if it exists. The LEP 
collaborations have all developed very similar searches for possible GMSB signals, as described 
in Refs. [^, |6^, The OPAL collaboration sets lower limits on the SUSY mass scale 

of A > 40,27,21,17,15TeV/c2 for N5 = 1,2,3,4,5 respectively, taking all NLSP lifetimes 
into account in a similar scan to the one described in this work. DELPHI quotes a limit of 
A > 17.5TeV/c^ for N5 = 4 and negligible NLSP lifetime. Only the ALEPH collaboration 
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has included neutral Higgs boson searches into the GMSB analysis, rendering better limits. 
As regards the NLSP mass, OPAL excludes them below 53.5, 87.4 and 93.7GeV/c^ in the 
neutralino, stau and slepton NLSP scenarios respectively. DELPHI has also looked for direct 
evidence of a massive sgoldstino S (the bosonic superpartner of the goldstino, with even Rp 
and produced via: e"'"e~ 6*7), in topologies with three energetic photons {S — 77) or two 
jets and one photon {S gg). This search allows direct exclusion of the SUSY breaking scale 
\/F^ up to 650 GeV for a light sgoldstino. 

To conclude. Fig. gives the status of the full MSSM spectrum in gauge mediated models 
after the LEP era. With the scan described in Chapter ^ and using only the ALEPH exclusions 
as described here, the final excluded masses at 95% confidence level have been computed for 
each messenger index N5. The ongoing LEP wide combination of results will improve these 
limits. 



6.3 Prospects 

The allowed parameter space in GMSB for searches where the stau is the NLSP and selectrons 
and smuons are lighter than the neutralino is extremely large, and will surely be tested in 
the near future with new data. At the incipient Run II of the Tevatron in Fermilab with 
= 2TeV, the dominant SUSY production is pp — > xtxi or X1X2 through off-shell W and 
Z. Charginos will tend to decay through their higgsino components by xf ~^ XiW^ and likewise 
for the neutralino X2 ~* Xi^*^. Thus for a neutralino NLSP scenario, final states consisting 
of W=^Z° + 77+ and W+W~ + 77+ will potentially lead to a discovery. The D0 
collaboration |^, expects to be able to discover charginos of masses up to 290 (340) GeV/c^ 
in the 77+ signature in the prompt x NLSP decays for a collected luminosity of 2fb~^ 
(30fb~^) by the end of the year 2004 (2007+). In the case of a stau NLSP, stau masses up to 
160 (200) GeV/c^ and chargino masses up to 340 (410) GeV/c^, will be discovered if the stau is 
long-lived for the same luminosities. The LEP limits {2a exclusion), as derived from the scan 
in this work, for the gluino and chargino masses, as a function of tan /3 and N5 can be seen in 



Fig. 3^. Overlaid in the chargino mass plot are the D0 discovery reaches (5ct) in the x NLSP 
case for short lifetimes and the stau NLSP case for long lifetimes. 

The task for the Large Hadron Collider at CERN (~2007) would then be to unravel the 
nature of SUSY breaking and measure the parameters of the theory with ~10% accuracy. 
Present and future colliders are opening a new energy range in physics, thus at the beginning 
of this century we are about to test new theories that may reshape our understanding of 
spacetime. . . or even better, be surprised by something completely unexpected. Whether it is 
supersymmetry or something else, it is surely going to be exciting. 
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Figure 6.1: Lower mass limit for all supersymmetric particles as a function of the number of messenger 
families. The limits are derived from the scan after applying GMSB and neutral Higgs boson searches 
results. Higher masses may be excluded for certain regions of parameter space. Thus these limits 
define the minimum allowed values for the masses for all possible GMSB models inside the ranges 
described in the scan (Tab. 5.1). 
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Figure 6.2: 95% C.L. exclusion on tlie gluino mass for (a) positive and (b) negative /i, for models with 
one or five messenger families, as derived from the GMSB searches. The lower limit on the chargino 
mass at 95% C.L. is displayed for (c) short and (d) long NLSP lifetimes respectively. The dotted lines 
represent the lower limit from GMSB and Higgs boson searches for N5 = 1 and the dashed line for 
N5 = 5. 
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